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Abstract 

We review and extend the theory of the dynamics of Bose-Einstein condensa- 
tion in weakly interacting atomic gases. We present in a unified way both the 
semiclassical theory as well as the full quantum theory. This is achieved by 
deriving a Fokker-Planck equation that incorporates both the coherent and 
incoherent effects of the interactions in a dilute Bose gas. In first instance we 
focus our attention on the nonequilibrium dynamics of a homogeneous Bose 
gas with a positive interatomic scattering length. After that we discuss how 
our results can be generalized to the inhomogeneous situation that exists in 
the present experiments with magnetically trapped alkali gases, and how we 
can deal with a negative interatomic scattering length in that case as well. 
We also show how to arrive at a discription of the collective modes of the 
gas that obeys the Kohn theorem at all temperatures. The theory is based 
on the many-body T- matrix approximation throughout, since this approxi- 
mation has the correct physical behavior near the critical temperature and 
also treats the coherent and incoherent processes taking place in the gas on 
an equal footing. 

PACS numbers: 03.75.Fi, 32.80.Pj, 42.50.Vk, 67.65.+Z 
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I. INTRODUCTION 



In the second half of the 1920's the founding fathers of quantum mechanics were mainly 
concerned with the understanding of the physical world on an atomic scale. In particular 
the study of the spectra of atomic hydrogen and helium led, besides the formulation of the 
famous Schrodinger equation Jl| and the interpretation of the wavefunction as a probability 
amplitude g], to the exclusion principle || or more generally the (anti)symmetrization of the 
wavefunction in the case of particles with (half)integer spin M. This important postulate 
of modern quantum theory divides all particles into two distinct families, which are now of 
course known as fermions and bosons, respectively. 

Already before the discovery of the (anti)symmetrization requirement, it was shown on 
statistical grounds only that under extreme conditions ideal gases of these indistin- 
guishable particles have also remarkable features on a macroscopic scale. Furthermore, in 
the case of an interacting gas even more stunning phenomena may occur, however weak the 
interaction. A well-known example in this context is superconductivity. In a superconduct- 
ing metal an attractive interaction between two electrons with opposite momenta causes an 
instability of the Fermi-surface and the formation of Cooper pairs || . The latter are allowed 
to move freely through the lattice resulting in a superconducting current and a vanishing 
resistance. According to the sucessful BCS-theory || describing this, the attractive inter- 
action is the result of a phonon exchange process. It is, however, believed that due to the 
weakness of this interaction only critical temperatures up to about 30 K are obtainable. 
Therefore, the underlying mechanism for the high-temperature superconductors cannot be 
explained in this way and is today still a very active field of investigation. 

Although the BCS-theory shows that superconductivity can in a certain sense be regarded 
as a result of a Bose-Einstein condensation |]6| JT0| , pTl] of the Cooper pairs, the most striking 
example of this condensation process was until recently always associated with the superfluid 
phase of liquid 4 He |12|]. Phenomenologically, the characteristics of superfluidity can be 



explained if the dispersion relation of the elementary excitations differs from the particle- 
like dispersion e(k) = h 2 k 2 /2m and is linear for small momenta [13] . To see that this feature 



indeed leads to superfluidity, we consider the situation of liquid He in a very long cylindrical 
pipe moving with velocity v along its symmetry axis Jl4]] . Describing the strongly interacting 



Bose liquid in the pipe as an assembly of noninteracting quasiparticles in a state of thermal 
equilibrium, we find from the usual methods of statistical mechanics that in the laboratory 
frame the number of quasiparticles with momentum Tik and energy huj(h) is given by 

^(k) = e (ft u ,(k)-ftk-v)/fe B T _ l ' (1) 

The velocity v is still arbitrary at this point, but its magnitude has an upper bound 
because the occupation numbers iV(k) must be positive. Therefore, we require for all k ^ 
that hu(k) > Kk. ■ v. In this case the total momentum P = J2k frkiV(k) carried along with 
the walls of the pipe, is at low temperatures T clearly much smaller than the momentum 
Nmv that we obtain when the whole liquid containing N atoms is moving rigidly with 
the walls. Hence, if the dispersion relation is linear, i.e., huj{k) = hck, we conclude by a 
Galilean transformation that the fluid can have a stationary (frictionless) flow if the velocity 



is small enough and obeys v < c Herewith, both the existence of superfluidity as well 
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as a critical velocity above which the phenomenon cannot take place is explained. It is 
important to note, however, that the above argument is semiclassical and predicts that a 
superflow is absolutely stable if the condition v < c is satisfied. This is incorrect in principle, 
because the superflow can actually decay by means of quantum and thermal fluctuations and 
is therefore only metastable |16[]. As we will see below, the effect of quantum and thermal 
fluctuations on the semiclassical picture is a recurrent theme throughout the following. 

At zero temperature, the required linear behavior of the dispersion relation was first 
derived on the basis of a microscopic theory by Bogoliubov [[17]] and was subsequently put 
on a firm basis by Beliaev [ffH], Hugenholtz and Pines [O], and Gavoret and Nozieres 



Simultaneously, it was also shown to be valid at nonzero temperatures by Lee and Yang [21 



using the pseudopotential method [p2]| , and more rigorously by Hohenberg and Martin |23 



and Popov [^4]] using field-theoretical methods. At this point it is important to mention also 
the work of Nepomnyashchh and Nepomnyashchii, who explicitly showed that the linear 
dispersion of the elementary excitations also survives if the infrared divergencies in the 
perturbation theory are taken into account exactly |25 . 



From these famous papers the linear dispersion at long wavelengths in a Bose system 
with short-range interactions is therefore well established. However, to actually calculate its 
velocity c from first principles is not feasible in general and in particular not for the strongly 
interacting 4 He liquid. To achieve that we have to consider a dilute Bose gas, for which we 
can rely on the weakness of the interactions or more precisely on the smallness of the gas 
parameter (na 3 ) 1 / 2 , where n is the density and a the positive s-wave scattering length of 
the particles in the gas. Unfortunately however, until recently Bose-Einstein condensation 
in a gas was never realized experimentally and the theory of the weakly interacting Bose 
gas [|D],[n|^l],|4[] could therefore not be applied to any real system. This situation changed 
dramatically in the summer of 1995 when within a few months three experimental groups 
reported the observation of Bose-Einstein condensation in a vapor of spin-polarized 87 Rb 
7 Li J27|, and 23 Na |2§ atoms. 
Historically, the first step towards the achievement of Bose-Einstein condensation in 



a dilute gas was made in 1980 by Silvera and Walraven ||29|| , when they demonstrated 
the stabilization of spin-polarized atomic hydrogen. After this pioneering experiment, the 
combined efforts of several experimental and theoretical groups around the world led to a 
detailed understanding of the atomic and statistical properties of this new quantum system 
| 30| . Moreover, just before the reports on the above mentioned experiment with the rubidium 
vapor, doubly spin-polarized atomic hydrogen gases still held the record in achieving the 
required high densities and low temperatures and were generally considered to be the most 
promising candidates for the observation of Bose-Einstein condensation [RTj]. Although, 
the latter expectation has clearly been proven wrong, it is fair to say that the sucessful 
experiments with the alkali vapors would not have been possible without the knowledge 
acquired by the atomic hydrogen community. In particular the development of magnetic 



traps [32j33| and evaporative cooling techniques p4|j35H has turned out to be of crucial 



importance for the actual achievement of Bose-Einstein condensation [STS] 



At present, the most important reason for the interest in degenerate Bose gases is that 
we can study the nonequilibrium dynamics of these trapped gases in such detail that experi- 
mental results can be directly confronted to microscopic many-body theories. In particular, 
the dynamics of condensate formation and the collective modes of an inhomogeneous Bose 
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condensed gas are two topics that have received considerable attention in the last three years 
and on which the present paper is also focussed. More generally, however, the aim of this 
paper is to develop with the techniques of quantum field theory a unified nonequilibrium 
theory for trapped atomic gases that is capable of treating the dynamics of the gas, not only 
in the normal and superfluid phases, but also during the phase transition itself. In agreement 
with the above goal of the paper, we refrain from reviewing other approaches that can and 
actually have been applied to the topics of interest here. Nevertheless, we have included a 
rather large amount of references. In this manner we aim to give, to the best of our knowl- 
edge, proper reference to important contributions to this field of physics. Furthermore, we 
also want to provide the reader with additional (and sometimes more advanced) literature 
that may be useful for a deeper understanding of various technical details, even though this 
is not necessary for a first reading, since the paper is intended to be self contained and 
essentially only assumes a working knowledge of the methods of second quantization, the 
path-integral formulation of quantum mechanics and the Hartree-Fock theory of interacting 
many-particle systems. 

The theoretical challenge that is posed by trapped atomic Bose gases can be stated as 
follows. On the one hand the interactions between the atoms influence the one-particle 
wavefunctions and energies of the gas. On the other hand the same interactions also change 
the occupation numbers of these states. As a result of the strong coupling between these, 
respectively, coherent and incoherent processes, the nonequilibrium dynamics of the gas is 
in principle quite complicated and requires a quantum kinetic theory that simultaneously 
incorporates both effects of the interactions. Such a quantum kinetic theory has recently 
been developed by Gardiner and Zoller on the basis of a perturbative treatment of the 
master equation for the many-body density matrix |[37]. Below we ultimately also arrive at 
a quantum kinetic theory, but this is achieved by using field-theoretical tools to derive a 
single Fokker-Planck equation that nevertheless describes both the coherent and incoherent 
dynamics of the gas. In principle we could, of course, also first derive a master equation 
and subsequently the same Fokker-Planck equation by making use of well-known procedures 
from quantum optics. However, we believe that the approach presented here is more direct 
and, most importantly, much more convenient if we systematically want to go beyond the 
pseudopotential method that is not sufficiently accurate for our purposes and that is used 
by Gardiner and Zoller. 

The paper is organized as follows. In Sec. |IJ we start with a discussion of an ideal Bose 
gas in a trap. In Sec. [11 A| we consider the isolated case, but in Sec. |11 B| we couple the gas 
in the trap also to a thermal reservoir. The reason for treating the ideal Bose gas is that in 
this way we are able to introduce most easily all the mathematical tools that are required 
for the discussion of an interacting Bose gas. The necessary extensions for an interacting 
gas are then developed in Sees. HJ, lYl an d 0- in these sections we mainly focus on the 
problem of understanding the dynamics of condensate formation, because the theory needs 
to be the most accurate and therefore the most general for this problem. In Sees. |T| and 
|TV| we first treat a homogeneous Bose gas and in Sec. [V] we then turn to the inhomogeneous 
case. Discussing the homogeneous Bose gas first is a useful intermediate step, because as 
a result of the translational symmetry the coherent dynamics of the gas involves only the 
one-particle energies and not their wavefunctions. Nevertheless, the dynamics is nontrivial. 
To get a better understanding of the timescales involved in Bose-Einstein condensation we 
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therefore first consider in Sec. [TTT| the semiclassical dynamics of the gas, in which fluctuations 
of the order parameter are neglected. Roughly speaking, the semiclassical theory is the 
nonequilibrium generalization of the Bogoliubov theory. Next, we discuss in Sec. [IV] the full 
quantum theory that incorporates also fluctuations of the order parameter. 

Sec. [V] is devoted to the quantum kinetic theory for inhomogeneous Bose gases. After 
our experience with the homogeneous Bose gas we can easily derive the desired Fokker- 
Planck equation for an interacting Bose gas trapped in an external potential. In contrast to 
homogeneous gases there are now, however, two cases to consider. If the density of the gas 
is sufficiently small, the average interaction between the atoms is always smaller than the 
energy splitting between the one-particle states in the trap and the coherent dynamics of the 
gas can be neglected. This most simple situation corresponds to the weak-coupling limit and 
is presented in Sec. |V A| . Experimentally, the weak-coupling limit is essentially only realized 
for atomic 7 Li due to the effectively attractive interatomic interactions of this element. As 



an application of the weak-coupling theory, we therefore discuss in Sec. |V B| the dynamics of 



Bose-Einstein condensation in a gas with effectively attractive interations in some detail. It 
is only after that, that we turn our attention to the most complicated strong-coupling limit 
and discuss three other applications of our general nonequilibrium approach. Since it follows 
most naturally from the theory presented in Sec. |V B| , we first briefly consider in Sec. |V C 



an interesting (zero temperature) property of a condensate, namely the so-called 'diffusion' 
of its phase. In Sec. [VD| we then return to the dynamics of Bose-Einstein condensation and 



in Sec. [VE| we finally present also a theory for the collective modes of a trapped Bose gas. 



We end the paper in Sec. [VI| with some conclusions and an outlook. 

II. CALDEIRA-LEGGETT LIKE TOY MODEL 

Before we start with the full problem of Bose-Einstein condensation in an interacting gas, 
we first consider in this section the case of an ideal Bose gas. In particular, we consider an 
isolated, noninteracting Bose gas in an external trapping potential V ex (x) with one-particle 
states Xo( x ) an d one-particle energies e a . In this manner we try to give an introduction 
to the Schwinger-Keldysh formalism for nonequilibrium processes [0,0 that we are going 
to use throughout this paper. A review of the Schwinger-Keldysh formalism has already 



been given by Danielewicz P0 |, but he uses the operator language which turns out to be 
less convenient for our purposes. Instead, we need a formulation in terms of functional 
integrals ||41||. In addition, we also couple this system to a reservoir in a way that is familiar 



from the Caldeira-Leggett model for a particle experiencing friction ||42|| . This gives us the 
opportunity to explain how dissipation, and as a result relaxation, is incorporated in the 
Schwinger-Keldysh formalism. This is especially important in Sec. |V], when we discuss the 
effect of thermal fluctuations on the semiclassical theory of Bose-Einstein condensation. 

A. Isolated Ideal Bose Gas 

In textbooks an ideal Bose gas is generally discussed in terms of the average occupation 



numbers of the one-particle states Xq( x ) [0. Given the density matrix p(t ) of the gas at 



an initial time t , these occupation numbers obey 
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N a (t)=Tr[p(t )ft a (t)4> a (t)\ , (2) 

with fojyt) and ip a (t) the usual (Heisenberg picture) creation and annihilation operators of 
second quantization respectively. Because the hamiltonian of the gas 

H = J2^ a (t)Mt) (3) 

a 

commutes with the number operators N a (t) = ■0^(t)'0 a ,(t), the nonequilibrium dynamics of 
the system is trivial and the average occupation numbers are at all times equal to their value 
at the initial time to- If we are also interested in fluctuations, it is convenient to introduce 
the eigenstates of the number operators, i.e., 

IW};*) = I P aK J> |o), (4) 

and to consider the full probability distribution 

P({N a };t) = Tr[p{t )\{N a y,t)({N a y,t\] , (5) 

which is again independent of time for an ideal Bose gas. The average occupation numbers 
are then determined by 

N a (t)= £ N a P({N a };t) (6) 

{N a } 

and the fluctuations can be obtained from similar expressions. 

As mentioned in the introduction, we are in the following also interested in the phase 
dynamics of the gas. To study that it is necessary to introduce so-called coherent states, 
which are eigenstates of the annihilation operators ip a (t) and are much used both in the field 
of quantum optics E| as well as condensed matter physics |46| . For our purposes it is even 



more convenient to consider eigenstates of the field operator 

^(X,t) =E^(*)Xa(x) . (7) 

a 

An eigenstate of V>(x, t) with eigenvalue 0(x) = Z) a 0aX«( x ) is given by 

\frt) = exp || dx 0(x)^( X) t)| |0> = exp j$>^!w} |0> (8) 

and is clearly also an eigenstate of ip a (t) with eigenvalue 4> a . The latter can also be proved 
immediately by noting that the equal-time commutation relation [ip a (t), ip^ (t)} = 1 implies 
that ip a {t) acts as d / dip^t) on these states. Moreover, from their definition it can easily be 
shown that these eigenstates obey the inner product 



exp {I rfx 0*(x)0'(x)} = exp 0>x} 



(9) 
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and the completeness relation 

\M)(M\ 



(M\M) 

if the integration measure of the functional integral is taken to be 



(10) 



In analogy with the occupation number representation in Eq. (|5|), we can now obtain a 
completely different description of the Bose gas by making use of the above coherent states 
and considering the probability distribution 



P[0*,0;t]=Tr 



\ ( .,\M)(M\ 



(M\M) 



(12) 



Although we expect that this probability distribution is again independent of time, let us 
nevertheless proceed to derive its equation of motion in a way that can be generalized in 
Sees. [Til] and [TV] when we consider an interacting Bose gas. First, we need to expand the 
density matrix p(t ) in the above coherent states. For an isolated Bose gas it is appropriate to 
take an initial density matrix that commutes with the total number operator N = J2 a N a (t) 
and we then find the expression 



2 , i \4>o;t )((j)o;t \ 



p(t ) = / «M0o] p^m-m , ' ' : \; , (is) 

in which the expansion coefficients p[|0o| 2 ;^o] only depend on the amplitude of the field 
4> (x) but not on its phase. This is equivalent to saying that the initial state of the gas does 
not have a spontaneously broken U(l) symmetry. Since we are ultimately interested in the 
dynamics of Bose-Einstein condensation it is of course crucial not to consider an initial state 
in which this symmetry is already broken and we therefore never use such initial states in 
the following. 

Next, we substitute this expansion into Eq. (|T2|) to obtain 

This is a particularly useful result, because the time dependence is now completely de- 
termined by the matrix element ((f); t\<f>o] to) for which the functional integral representa- 
tion is well known [f46fl . It is given by a 'path' integral over all complex field evolutions 
if)(x.,t+) = J2 a '4 ) a{t+)Xa{' x ) from to to t. More precisely we have 



t\fat ) = f (x,t) * {x) d[r}dm exp{^s + [rM} > (1 5 ) 

Jib(x,t )=<t>o(x.) in ) 



rt/.*(x,t)=r(x) ^ n ( i 

/i/)(x,t )=<Ao(x) 

with the forward action S + [ip*, tp\ given by 

«■* . / . 



S + [r, V] = E { -iWaWaW + f dt+ F a {t+) - e a \ MU) I • (16) 
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In the same manner the matrix element (0; t|</>o; to)* = (0oj to\4>', t) can be written as a 'path' 
integral over all field configurations ^(x, t_) = J2a V ; a(^-)Xa( x ) evolving backward in time 
from t to t , i- e -> 

<0;t|0 o ;to>* = / exp \-S-[rM \ , (17) 

with a backward action 

= E l-iW a (t)Mt) + jf° dt- Mt-) - e «) • (is) 



Putting all these results together, we see that the probability distribution P[<j>*, (ft; t] can 
in fact be represented by a functional integral over all fields ^(x, t) that evolve backwards 
from t to to an d then forward in time from to to t. Absorbing the factor p[|0 o | 2 ; to] into the 
measure of the functional integral, we thus arrive at the desired result that 

/ -V*(x i *)=0*(x) f I -I 

PW, 0; t] = / ^ d[r\dm exp -S[P, ^ , (19) 

Ji/>(x,t)=0(x) l n ) 

where the total (backward-forward) action in first instance obeys 

s[p,i>] = S-WM + s+[rM = -^E (C(t)^(t) - |</g 2 ) 

a 

+T,J ct dt ' \^^')^Mt') - , (20) 

and the integration along the Schwinger-Keldysh contour C l is defined by J* c < dt' = j/° dt_ + 
/ t * <it + p8|^0| . Note also that in Eq. (|19|) we have explicitly specified the boundary condi- 
tions on the functional integral. It is interesting to mention that these boundary conditions 
are essentially dictated by the topological terms in the action S[i/)*,ip], which is a gen- 
eral feature of the path-integral formulation of quantum mechanics due to the fact that 
the quantum theory should have the correct (semi) classical limit Making use of the 

periodicity of the field i[)(x.,t) on the Schwinger-Keldysh contour, the variational principle 
5S[ip*, ip]/8ipa (t±) = indeed leads not only to the Euler-Lagrange equation 

d 

ih-^-ip a (t ± ) = e a ?p a (t±) , (21) 

which agrees with the Heisenberg equation of motion ihdip a (t±) / dt± = [tf) a (t±),H] for the 
annihilation operators, but also to the appropriate boundary condition Sip^(t) = 0. In the 
same way we find the complex conjugate results if we require that 5S[i/;*, if)]/6if) a (t±) = 0. 
Substituting the boundary conditions in Eq. ([20) and performing a partial integration, we 
then finally obtain for the action 

S[P, = £ j ct dt' C(t') - e a ) Mt') , (22) 
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which we for completeness sake also rewrite as 



S\rM = J ct dt'J dM*(x,0 + v*(x)^(x,0 • (23) 

We are now in a position to derive the equation of motion, i.e., the Fokker-Planck equa- 
tion, for the probability distibution P[0*,0;i]. This is most easily achieved by performing 
the variable transformation ip(?L,t±) = 0(x, £') ± £(x, t')/2 in Eq. (|T9|). In this manner the 
fields ?/>(x, £_) and ^(x, t + ) that live on the backward and forward branch of the Schwinger- 
Keldysh contour, respectively, are 'projected' onto the real time axis. Moreover, at the same 
time we effect a separation between the (semi) classical dynamics described by 0(x, t') and 
the quantum fluctuations determined by £(x, t'). After the transformation we have 

PWA;t] = f* (X ' )=0 (X) rf[0*]d[0] [ d[CW\ exp{~S[«f ,</>;£*,£]} , (24) 
J <f>(x.,t)=<t>(x) J L tl ) 

with 

^[0*,0;r,e] = E//^ {c(o(^|7- e ^e Q (O+e:(o(^|7- e ^0a(o} • (25) 

Because this action is linear in £ a (t') and the integration over the quantum fluctuations 

leads just to a constraint and we find that 



d e n \ .*,., s ( . d e 



(26) 



P[0*,0;t] = / d[<P*}d[<P} US 

J0(x,i)=(/)(x) „ 

or equivalently that 

P[0*,0;t] = / c#SM0o] P[|0o| 2 ;io]II 5 (l^-^)| 2 ) ' ( 2? ) 

J a 

where we used that P [</>*, 0; to] is only a function of the amplitude |0| 2 and also introduced 
the quantity 0„(t) obeying the semiclassical equation of motion 

ih^tit) = e a <Pt(t) (28) 

and the initial condition 0^(to) — 4>o,a- 

The latter equation is thus solved by 0„(t) = (f) 0:a e~ lta ^~ to ^ h and we conclude from a 
simple change of variables in Eq. ( pTf) that for an ideal Bose gas P [</>*, = P[|0| 2 ; to], as 
expected. We also see from Eq. ( PT| ) that the desired equation of motion for P [</>*, <p; t] reads 

ih^P[<P, M = ~ (E ^ e ^) & *1 + (E ^T e ^) ^ *1 • ( 29 ) 

This is indeed the correct Fokker-Planck equation for an ideal Bose gas ||8| . To see that we 
first consider the average (<j> a )(t) = J d[<p*]d[<p] <p a P[(f)*,<p;t]. Multiplying Eq. (|^) with <p a 
and integrating over 0(x), we easily find after a partial integration that 
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i^{Mt) = e a (<f> a )(t) , (30) 

which precisely corresponds to the equation of motion of (i(j a (t)) = Tr[p(t Q )ilj a (t)} in the 
operator formalism. Similarly, we find that 

ft^ifaKt) = -eaifcKt) , (31) 

in agreement with the result for (4>a{t)) = Tr[p(t )^(t)]. 

Next we consider the average of \4> a \ 2 , for which we immediately obtain 

ih^{\<f> a \ 2 )(t) = 0. (32) 

We expect this result to be related to the fact that in the operator formalism the occupation 
numbers N a (t) are independent of time. Although this turns out to be true, to give the 
precise relation between (\4> a \ 2 ){t) and N a (t) is complicated by the fact that at equal times 
the operators ip a {t) and ipa(t) do not commute. However, the path-integral formulation of 
quantum mechanics is only capable of calculating time-ordered operator products [ |49(| . In 
our case this implies that (\4> a \ 2 ){t) is the value at t' — t of 



Tr 



)(t )T c t (4(*)&(0)] = ®(t,t')Tr [p(f )4(*)&(0] + ®(t',t)Tr [p(* )$,(04(*) 



with Tct the time-ordening operator on the Schwinger-Keldysh contour and Q(t,t') the 
corresponding Heaviside function. Since the Heaviside function is equal to 1/2 at equal 
times we conclude that 

^ a \*)(t) = N a {t) + l - (33) 

and that Eq. (|32| ) is thus fully consistent with the operator formalism. An intuitive under- 
standing of the relation between (\cf) a \ 2 )(t) and the occupation numbers can be obtained by 
noting that, since all our manipulations up to now have been exact, we expect (|0 Q | 2 )(t) to 
contain both classical and quantum fluctuations. These correspond precisely to the contri- 
butions N a (t) and 1/2, respectively. 

Finally we need to discuss the stationary solutions of the Fokker-Planck equation. It 
is not difficult to show that any functional that only depends on the amplitudes \<f) a \ 2 is 
a solution. As it stands the Fokker-Planck equation, therefore, does not lead to a unique 
equilibrium distribution. This is not surprising, because for an isolated, ideal Bose gas there 
is no mechanism for redistributing the particles over the various energy levels and thus for 
relaxation towards equilibrium. However, the situation changes when we allow the bosons in 
the trap to tunnel back and forth to a reservoir at a temperature T. The corrections to the 
Fokker-Planck equation that are required to describe the physics in this case are considered 
next. However, to determine these corrections in the most convenient way, we have to 
slightly generalize the above theory because with the probability distibution P [</>*, </>;£] we 
are only able to study spatial, but not temporal correlations in the Bose gas. 



To study also those we follow a well-known procedure in quantum field theory pi. 52 



and construct a generating functional Z[J, J*} for all (time-ordered) correlation functions. 
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It is obtained by performing two steps. First, we introduce the probability distribution 
Pj[4>* , 0; t] for a Bose gas in the presence of the external currents J(x, t) and J*(x, t) by 
adding to the hamiltonian the terms 

-nfdx ($(x,f)j*(x,t) + j(x,t)^t( x ,t)) = -hJ2(Mt)J* a (t) + J a ml(t 

As a result we have 

rV*(x,t)=0*(x) 



Pj[0*, 0; t] = / d[rW\ exp -S^*, ^] 

Jip(x,t)=4>(x) in 



exp^ij dif J dx (^(x,t')^*(x,t') + J(x,t')V'*(x,t / )) 



(34) 



Second, we integrate this expression over 0(x) to obtain the desired generating functional. 
Hence 



Z[J,J*] = / 



xexpjz J dt J dx (^(x, t) J*(x, t) + J(x, t)4>*(x, t)) 



(35) 



It is important to realize that Z[J, J*] is indeed independent of the time t because of the 
fact that Pj[(j)*,(j)]t] is a probability distribution (cf. Eqs. ( |10|) and dl2|) ) and thus properly 
normalized. We are therefore allowed to deform the contour C* to any closed contour that 
runs through to- Since we are in principle interested in all times t > to, the most convenient 
choice is the countour that runs backward from infinity to to and then forwards from to to 
infinity. This contour is denoted by C°° and also called the Schwinger-Keldysh contour in 
the following because there is in practice never confusion with the more restricted contour 
C l that is required when we consider a probability distribution. With this choice it is also 
clear that Z[J, J*} is a generating functional. Indeed, Eq. (|35j) shows explicitly that all time- 
ordered correlation functions can be obtained by functional differentiation with respect to 
the currents J(x, t) and J*(x, t). We have, for instance, that 



(36) 



J,J*=0 



and similarly that 



Tr 



p(to)T C cc(^(x,t)^t( x ',t') 



8 2 



i 2 SJ*(x,t)5J(x.',t 



(37) 



j,j*=o 



Note that the times t and t' always have to be larger or equal to to for these identities to be 
valid. 



B. Ideal Bose Gas coupled to a Reservoir 

For the reservoir we take an ideal gas of N bosons in a box with volume V. The states 
in this box are labeled by the momentum fok and equal to Xk( x ) = e* kx /v^- They have 
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an energy e(k) = h 2 k 2 /2m + AV ex , where AV ex accounts for a possible bias between the 
potential energies of a particle in the center of the trap and a particle in the reservoir. The 
reservoir is also taken to be sufficiently large that it can be treated in the thermodynamic 
limit and is in an equilibrium with temperature T and chemical potential /i for times t < to. 
At to it is brought into contact with the trap by means of a tunnel hamiltonian 

H mt = 4tEE (ak)4wto + c(k)&(*)&(*)) . (38) 



V 

with complex tunneling matrix elements t a (k) that for simplicity are assumed to be almost 
constant for momenta Hk smaller that a cutoff hk c but to vanish rapidly for momenta larger 
than this ultraviolet cutoff. Moreover, we consider here only the low-temperature regime 
in which the thermal de Broglie wavelength A = (2-kTi 2 /inksT) 1 / 2 of the particles obeys 
A ^> l/k c , since this is the most appropriate limit for realistic atomic gases. 

To study the evolution of the combined system for times t > to we thus have to deal 
with the action 

s[r,^r R ,M = -4tEE / dt (t a (k)vv,(t)V£(*) +t* a (k)Mt)r a (t)) 

V V a k C °° 

+ E j CBB dt C W -e a + f)j Mt) + E Xoc dt ( in M - e(k) + ^ ' (39) 

if we measure all energies relative to the chemical potential and also introduce the complex 
field t/)_r(x, t) = J2k '0k(^)Xk( x ) for the degrees of freedom of the reservoir. However, we 
are only interested in the evolution of the Bose gas in the trap and therefore only in the 
time-ordered correlation functions of this part of the system. The corresponding generating 
functional 

z[j,r] = J d[rw\ J d[r R w R \ ex P {^*,^«,^]} 

xexpji^ dtjdx (^(x,t)J*(x,t) + J(x,t)^*(x,t))| . (40) 

is of the same form as the functional integral in Eq. (|35|), but now with an effective action 
that is defined by 

exp j^ e// [^*,^]} = J d[r R w R \ exjp{^s[r,ip;r R ,ipR)} . (41) 

Hence, our next task is to integrate out the field t/>_r(x, t), which can be done exactly because 
it only requires the integration of a gaussian. 

To familiarize ourselves with the Schwinger-Keldysh formalism, we perform the gaussian 
integration here in some detail. In principle, this can of course be done explicitly by making 
use of the fact that for the initial density matrix p R (to) of the reservoir we have that 

Pn[\M 2 M-Jlj^f Mmii) , (42) 

with iV(k) = l/(e^( e w - w _ i) the appropriate Bose distribution function and (3 = 1/ksT. 
However, in practice it is much more easy to use a different procedure. It is based on the 
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observation that if we introduce the 5 function on the Schwinger-Keldysh contour defined 
by/c°° dt' ^(t,t') = 1 and the Green's function G(k;t,t') obeying 



(ih^ - e(k) + f?j G(k; t, t') = M(t, t') , 



(43) 



the action S[ip*, ip; tjj R , i^r] can be written as a complete square, or more precisely as the sum 
of two squares Si[ip* , ip] and S2I1P*, ip; ipRi 4>r] that are given by 

-TT, EE/ dtj dt' r a (t)t*MG(k-,t,t')t a ik)^(t') (44) 



a, a' k 

and 



S 2 [r, ^ p R , M = E / coo * (ft® - J^y E \ Qoa dt' t*Mr*(t')G(k; t', t) 



4=E/ dt' G{^t,t')^ a {t')t a {^)\ , (45) 



x - e(k) + f?j L w (t) - - 

respectively. Since the first term is independent of the field iPr(x., t), we only need to 
evaluate the functional integral J d[ip R ]d[ip R ] exp (iS 2 [ip*, ip; ip R , ip R ]/h). Performing a shift 
in the integration variables, we however see that this functional integral is equal to 

j d[r R W R ] expj^E^^^kW U^-e(k)+^ k (t)\=Tr[p R (t )] = l . (46) 

As a result the effective action S e ^[ip*, ip] is just equal to Si[ip*, which can be slightly 
rewritten to read 

S eff [ifj*,*fj} 

= E J coo dtj c<x dt' r a (t) I (ih^ -ea + Jj 5 a , a ,8(t 7 t') - fr£ a , a ,(M')} <M0 , (47) 

with the so-called selfenergy E Q Q ,/ (t, t') obeying 

hX aja ,{t,t') = ^XX(k)G(k;M%'(k) ■ (48) 

This is our first example of an effective action describing the dynamics of a Bose gas. 
Before we can study its consequences we clearly first need to determine the Green's function 
G(k; t, t') in terms of which the selfenergy is expressed. Although we know that this Green's 
function fulfills Eq. (fE|), we cannot directly solve this equation because we do not know 
the appropriate boundary condition at t = t'. To calculate G(k;t,t') we therefore have to 
proceed differently. It is however clear from Eq. ([13|) that G(\t;t,t') is a property of the 
reservoir and we thus expect that it can somehow be related to a time-ordered correlation 
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function of this reservior. To see that explicitly we consider again the generating functional 
of these correlation functions, i.e., 

Z R [J,J*]= Jd[r R ]d[^ R ] expjl^^d* Vk(*) ^-e(k)+/iUkW} 

xexp!^J2j coo dt (Mt)JZ(t) + Mt)i>l(t))j ■ (49) 

It is again a gaussian integral and can thus be evaluated in the same manner as before. The 
result is now 



Z R [J, J*} = exp(-i£ j dt f dt 1 JZ(t)GQt;t,t')Mt') 

I k C °° C °° 



(50) 



which shows by means of Eq. (]37]) that 

iG(k;t,t') = Tr \p R (t )T c ~ U^t)^ 1 ) 



Note first of all that this is indeed consistent with Eq. 
obeys 



(51) 

because the right-hand side 



%h— Tr 

dt 



p R {ta)Tc°° [Mt)i>Ut' 



+ Tr 



ih5(t,t') + (e(k) - ti)Tr\p R (t )Tc 



( d ~ 



(52) 



in the operator formalism. Moreover, from this identification we see that the desired solution 
fulfilling the appropriate boundary conditions is apparently 



G(k; t, t') = _ ie -*(e(k)-*0(t-f )/» {e(t, f) (1 + N(k)) + G{t', t)N(k)} . 



(53) 



The specific dependence on the backward and forward branches of the Schwinger-Keldysh 
contour is thus solely determined by the Heaviside function Q(t,t'). As a result it is conve- 
nient to decompose the Green's function into its analytic pieces G > (k; t — t') and G < (k; t — t') 
by means of 



G(k; t, f) = 8(t, t , )C > (k; t-t') + 0(t', t)G < (k; t - t') . 



(54) 



Due to the fact that we are always dealing with time-ordered correlation functions, such 
a decomposition turns out to be a generic feature of all the functions on the Schwinger- 



Keldysh contour that we will encounter in the following ||54|| . For a general function F(t, t'), 
it is however also possible to have 5-function singularities. If that happens the correct 
decomposition is 



F(t, t') = F s (t)5(t, t') + 6(t, t')F > (t, t') + Q(t', t)F < (t, f) 



(55) 



This more general decomposition is not required here, but will be needed in Sees. |TT| and 
1V| when we determine the effective interaction between two atoms in a gas. 
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Having obtained the Green's function of the reservoir, we can now return to our dis- 
cussion of the effective action S e ^[if>*, if)] for the Bose gas in the trap. Although we have 
chosen to derive in Eq. ( f4?|) the effective action for the generating functional Z[J, J*], it is 
straightforward to show that the effective action for the probability distribution P[<f>* , (f>; t] 
is obtained by only replacing the contour C°° by C l . In the following we therefore no longer 
always specify all the boundary conditions on the time integration, if the precise details 
of this integration are not important and the discussion applies equally well to both cases. 
Keeping this in mind, we now again perform the transformation if> a (t±) = 4> a (t) ± £ a (t)/2 
to explicitly separate the (semi) classical dynamics from the effect of fluctuations. It leads 
in first instance to [E5| 



S e "[<f)\ 0; £*, £] = £ jf dt {^(t) (ihj t -e a + t?j at) + C(t) (ih^ t -e a + f?j Q (t)| 

-Y, f dt f df te(t)nsS^(t- 0^(0 + 0^(0) 

* Zj J to Jto v J 

- \ E / dt f dt' Cit^lAt - 0M0 , (56) 

act' 

where we introduced the retarded and advanced components of the selfenergy 

- = ±e> (±(* - 0) {KAt - - KAt - 0) (57) 

that affect the terms in the action that are linear in £*(£) and £ a (t), respectively, and also 
the Keldysh component 

KAt - = {KAt - + s a,^ - 0) (58) 

that is associated with the part quadratic in the fluctuations. 

The physical content of these various components of the selfenergy is understood most 
clearly if we now apply the beautiful functional-integral procedure that is due to Hubbard 
and Stratonovich [^6|,^7],[46|] . The basic idea is to write the factor 

ex p|-^E/ dtj dt' c a mlAt-t')tAt ! ) 

in the integrant of the functional integral / d[(p*]d[4>] J exp (iS eff [</)*, (f>; £*, t]/fi) 

as a gaussian integral over a complex field r/(x, t). It is equivalent, but in practice more 
convenient, to just multiply the integrant by a factor 1 that is written as a gaussian integral 
Jd[r]*]d[r]] exp(tS eff [r]%7]]/h) with 

s B/ V^] = ^E/ dtj dt' (2<(t)-W dt" C"(t")n^At"-t)) 

x (hZ K )-^(t - f 21^(0 - E / dt" hE« >a „(t - f)U(0) (59) 

V a" Jto I 

a complete square. Adding this to S e ^[4>*, 0; £*, £] the total effective action becomes 
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= E / dt J dt ' \ [ lh §- t + V*a(t)j $a,a>5(t - f) - KE^t - t')\ ^(f) 

+ E / d * / dt' £(t) | (fc - e a + n - jj a (t) j 5 a ^5{t - t') - KE^t - t')\ 



0a' (0 



+ 2 £ / <ft / eft' ifiCfXfts*)^ - 0*^(0 (60) 

a, a' 

and is thus linear in and £ Q (£). Integrating over these fluctuations we conclude from 
this action that the field 0(x, t) is constraint to obey the Langevin equations [58 

i^<i>a(t) = (e« - A*)0a(t) + E / ^ ^ S aa'(* ~ 0M0 + Va(t) (61) 
Ot i J t 

a' u 

and 

- i^lif) = (ea " /*)&(*) + E / ^ C^O^H^' " *) + »£(*) ( 62 ) 

with gaussian noise terms rj a (t) and 77* (t) that from the last term in the right-hand side of 
Eq. (|K)D are seen to have the time correlations 



(v:(t)vAt f )) = -fKoS - = \ J (^3 (! + 2iV ( k )) C(k)e-« k )^^')/X'(k) (63) 
in the thermodynamic limit. Moreover, in that limit the retarded selfenergy is equal to 

KE$,(t - = -8(t - i')| / ^3 C(k)e-^ k )^^')/^(k) (64) 

and the advanced selfenergy can be found from the relation E a , a (t — t) = (S^a' (t — t)j . 
Note that the Heaviside functions in these retarded and advanced selfenergies are precisely 
such that Eqs. (]6Tf) and d62|) are causal, i.e., the evolution of 0(x, £) and 0*(x, i) depends 
only on their value at previous times. This is clearly a sensible result and explains why these 
components of the full selfenergy enter into the terms in the action that are linear in the 
fluctuations. 

To understand why the Keldysh component enters into the terms that are quadratic in 
the fluctuations is more complicated and is related to the famous fluctuation-dissipation 
theorem [53,5^,5^]. Physically, the fluctuation-dissipation theorem ensures in our case that, 
due to the interactions with the reservoir, the gas in the trap now relaxes towards thermal 
equilibrium in the limit t — > 00. Let us therefore study this limit in somewhat more detail. 
However, before we do so we would like to mention for completeness that if we apply the 
above Hubbard- Stratonovich transformation to the functional integral for the probability 
distribution P[0*,0;t],we immediately find that analogous to Eq. 



P[<P,<j>;t] = J d[ V *]d[ V ] P[ V *,rj] J dttMfo] P[|0o|Vo]n^K-<^)| 2 ) . W 
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where now obeys the Langevin equation in Eq. (|6lD for a particular realization of 

the noise r) a (t) and with the initial condition 0„(to) — 4>o,a- Furthermore, the probability 
distribution for each realization of the noise is given by 

Ptf, V ] = exp || £ f dt> f dt" V* a (m^ K )a^(t' - t") Va ,{t") \ . (66) 

This result shows explicitly that, as expected, the Fokker-Planck equation associated with 
the Langevin equations in Eqs. fl6~Tf) and (|6~^) is in fact the desired Fokker-Planck equation for 
the probability distribution P [</>*, 0; £]. Hence, this Fokker-Planck equation can be obtained 
either from these Langevin equations or directly from the effective action in Eq. fl56|) , without 
performing a Hubbard-Stratonovich transformation. We now first follow the former route, 
but lateron also explain the more direct procedure because that is interesting in its own 
right and also shows the consistency of the theory. 

After this slight digression, we now return to the long-time behaviour of the gas. In 
the long-time limit two important simplifications occur. Looking at the expresions for the 
selfenergies, we see that their width as a function of the time difference t — t' is at most 
of order Ofi/hsT) in the low-temperature regime of interest where the thermal de Broglie 
wavelength of the particles obeys A ^> l/k c . If therefore t ^> to + 0(h/kBT) we are allowed 
to take the limit to ~~ > — oo. Taking this limit means physically that we neglect the initial 
transients that are due to the precise way in which the contact between the trap and the 
reservoir is made, and focus on the 'universal' dynamics that is independent of these details. 
In addition, at long times the dynamics of the gas is expected to be sufficiently slow that 
we can neglect the memory effects altogether. Finally, we consider here first the case of a 
reservoir that is so weakly coupled to the gas in the trap that we can treat the coupling with 
second order perturbation theory. As a result we can also neglect the nondiagonal elements 
of the selfenergies. In total we are then allowed to put (but see below) 

Z^«(t - if) ~ Sif(e tt - fi)5 a>a ,5(t - t') = ^> K 5 a , a ,5(t - t') , (67) 
where the Fourier transform of the selfenergies is defined by 

. „-ie(t-t')/% 

sS"(*-'') = /*ES K (<)^— ■ (68) 
Note that this implies in general that Yr a , a (e) = (S^, (e)j , and in particular therefore that 
With these simplifications our Langevin equations become 

ifr^<f>a(t) = (<* + ^ { a +) - + Va(t) (69) 

and just the complex conjugate equation for 0* (t). The retarded selfenergy in this equation 
is given by 
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with e+ the usual notation for the limiting procedure e a +iO. It clearly has real and imaginary- 
parts that we denote by S a and —R a , respectively. Denoting the Cauchy principle value 
part of an integral by V, they obey 

s - = J w*< k >^kj*-< k > (71) 

and 

r (7k 

R Q = ir J 5(e Q - e(k))\t a (k)\\ (72) 

The interpretation of these results is quite obvious if we consider the average of the Langevin 
equation, i.e., 

d 

ift-Qj;(<f>a)(t) = (e« + S a - iR a - fi)(<p a )(t) , (73) 

which is solved by 

((j) a )(t) = (0 Q )(O)e- i(e » +5 «-' i)t/R e- jRa * /?l . (74) 

Hence the real part of the retarded selfenergy S a represents the shift in the energy of state 
Xa( x ) due to the coupling with the reservoir. Indeed, Eq. ( ffll) agrees precisely with the 
value of this shift in second order perturbation theory. Moreover, the fact that \{4> a ){t)\ 2 = 
\{4> a )(0)\ 2 e~ 2Rat / h shows that the average rate of decay r a of the state Xo( x ) is equal to 
2R a /h, which, together with Eq. (72|), exactly reproduces Fermi's Golden Rule. 

From the equation of motion for the average (4> a ){t) we can also conclude that the right- 
hand side of the Fokker-Planck equation must contain the 'streaming' terms 



However, there must now also be a 'diffusion' term due to the fact that the average (\4> Q \ 2 ){t) 
is no longer independent of time since the interactions with the reservoir can lead to changes 
in the occupation numbers N a (t). To obtain the 'diffusion' term we thus need to determine 
ihd(\<j) a \ 2 )(t)/dt. To do so we first formally solve the Langevin equation by 

a (t) = e -i(e a+ h^-v)t/n |^ (0) _^J* dt , ^y^+ftEW-^f /ft j . ( 7 5) 

Multiplying this with the complex conjugate expression and taking the average, we obtain 
first of all 

(\<f> a \ 2 )(t) = e- 2R ^ h {(|0 Q | 2 )(O) + lx£ jf*df e^^'/ftj > (76) 
which shows that 

in^(m 2 ){t) = -2iR a {\<p a \ 2 )(t) - iftE? . (77) 
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Moreover, the Keldysh component of the selfenergy is given by 

hX« = -2m J ^ (1 + 2N(k))5(e a - e(k))|t a (k)| 2 (78) 
and therefore obeys 

KE* = -2i(l + 2N a )R a (79) 

with N a = (e^* 6 " - ^ — 1) _1 the Bose distribution function. This is in fact the fluctuation- 
dissipation theorem, because it relates the strength of the fluctuations determined by hT,^ , 
to the amount of dissipation that is given by R a . As mentioned previously, the fluctuation- 
dissipation theorem ensures that the gas relaxes to thermal equilibrium. This we can already 
see from Eq. fl77|), because substitution of the fluctuation-dissipation theorem leads to 

^<I<^| 2 >W = -2tR a (\<j) a \ 2 }(t) + iR a (l + 2N a ) , (80) 

which tells us that in equilibrium (|0 Q | 2 ) = N a + 1/2. In fact, we have argued in Sec. |11 A 
that in general (|0 Q | 2 )(t) = N a (t) + 1/2. Substituting this identity in Eq. (|30"D and realizing 
that T a N a = r a iV(k) due to the energy-conserving S function in T a , we indeed obtain the 
correct rate equation for the average occupation numbers 

^N a (t) = -T a N a (t) + T a N a = -T a N a (t)(l + N(k)) + T a (l + N a (t))N(k) , (81) 

that might justly be called the quantum Boltzmann equation for the gas, because the right- 
hand side contains precisely the rates for scattering into and out of the reservoir. 

Furthermore, Eq. (|77| ) shows that the 'diffusion' term in the Fokker-Planck equation is 

because the first term in the right-hand side of Eq. ( [77] ) is also present without the noise and 
is therefore already accounted for in the 'streaming' terms. In total we have thus arrived at 

iK^P[<i>*,frt] = 

~ (E t^(o* + - p)^ P[<f>\ 0; t] + (£ ^(e a + JffiH - ^ P[0*, 0; t] 

Using again the fluctuation-dissipation theorem, it is not difficult to show that the stationary 
solution of this Fokker-Planck equation is 

^*»]=n^^- P {-^-L^ w 4. (83) 
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Although this appears to be quite reasonable, it is important to note that this result is 
not fully consistent because the 'streaming' terms in Eq. fl8"2[) show that the energies of the 
states x«( x ) are shifted and equal to e a + S a . We therefore expect that the equilibrium 
occupation numbers are equal to the Bose distribution evaluated at those energies and not 
at the unshifted values e a . The reason for this inconsistency is, of course, that we have 
taken Sw ,Jf equal to EW ,,( (e a — /x), which was only justified on the basis of second-order 
perturbation theory. We now, however, see that to obtain a fully consistent theory we 
must take instead ^^' K = (e Q + S a — /j) . This implies that to obtain the retarded 

selfenergy, we first have to solve for the real part 

^/^^W^T^'-W (84) 

and then substitute this into 

r r/k 

Ra = 7T J J2^ 6{€a + S °~ e(k))Mk)|2 (85) 

to obtain the imaginary part. The fluctuation-dissipation theorem then again reads ftS^f = 
— 2i(l + 2N a )R a but now with the equilibrium distribution N a = ( e P( ea+Sa ~i J -) — l) -1 . As a 
consequence the quantum Boltzmann equation and the stationary state of the Fokker-Planck 



equation will now also contain these correct equilibrium occupation numbers [60 



This almost completes our discussion of this Caldeira-Leggett like toy model, that nev- 
ertheless already shows many of the features that we will encounter in Sec. |T| and |V] when 
we study a realistic interacting Bose gas. However, we promised to derive the Fokker-Planck 
equation also derictly from the effective action S e ^[<f)*, 0; £*, £], without making use of a 
Hubberd-Stratonovich transformation and the Langevin equations that are an immediate 
result of that transformation. Making the same approximations on the selfenergies as before, 
the effective action for the probability distribution P[<f>*,<j>;t] reads 

S e ^[0*, 0; T , £]=£ £ dt' C(t') (ih-^ -e a - ftE<f> + A*J Ut') 

+ E £ dt' £(t') U-jjp -e a - *£<+> + Q (t') -\Y,f to dt ' C(t')h^Ut') , (86) 

and is quadratic in the fluctuation field £(x, t). We can thus again perform the integration 
over this field exactly. The result is 



S eff W, 0] = E £ ^ [ih-jfi - t a - + /i j Ut') = £ dt' L(t') . (87) 



d 



rt 
Ho 



Since the probability distribution P[0*, 0; t] is equal to the functional integral 

P[0*, 0; t] = £ (X ' t)= * (X) d[0*M0] exp \^S ef f[P, 0] j , (88) 

we know from the general connection between classical mechanics and the functional for- 
mulation of quantum mechanics, that P[0*, 0; t] must obey the Schrodinger equation that 
results from quantizing the classical theory with the lagrangian L(t). 
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Fortunately, the quantization of this theory is straightforward. The momentum conjugate 

to (f> a (t) is 

*•<*> = m§m = m M - f - " SE " + ") m • m 

whereas the momentum conjugate to i.e., 7r*(t), is given by the complex conjugate 

expression. The corresponding hamiltonian is therefore 

= E ( e « + ^« +) - + i<V) ( e « + - A*) 

+ £^fk«(0l 2 - (so) 

Applying now the usual quantum mechanical recipe of demanding nonvanishing commuta- 
tion relations between the coordinates and their conjugate momenta, we can put in this case 
7Tq = (h/i)d/d(/) a and similarly 7r* = {%/i)d/d(j)* a . The Schrodinger equation 

ih^P[<f>*,4>;t} = HP[4>*,4>;t] (91) 



then indeed exactly reproduces the Fokker-Planck equation in Eq. fl82|). It should be noted 
at this point that the hamiltonian in Eq. Q9"H| ) has in principle ordening problems, due to the 
fact that the coordinates and the conjugate momenta no longer commute after quantization. 
These ordening problems are, however, resolved by realizing that in the path-integral formu- 
lation of quantum mechanics we always deal with so-called normally ordered hamiltonians 
in which the momenta are positioned to the left of the coordinates. For that reason we have 
to use the same ordening in Eq. (|90|). 

Having arrived at the same Fokker-Planck equation, we have now demonstrated the 
use of essentially all the theoretical tools that are required for a discussion of a weakly 
interacting Bose gas. Before we start with this discussion, however, we want to make a final 
remark about the effect of the nondiagonal elements of the selfenergies. Physically, including 
these nondiagonal elements accounts for the change in the wavefunctions Xa( x ) due to the 
interaction with the reservoir. It can clearly be neglected if the coupling with the reservoir 
is sufficiently weak, or more precisely if \tiE^ a (e a + S a — fi) \ is much smaller than the energy 
splitting \e a i — e a \ between the states in the trap. For a quantum dot in solid-state phyisics 
pl| this is often the case and for that reason we might call our toy model a bosonic quantum 
dot. However, in magnetically trapped atomic gases the average interaction energy between 
the atoms is already slightly below the Bose-Einstein condensation critical temperature 
larger than the typical energy splitting and the nondiagonal elements of the selfenergies 
are very important. Such a strong-coupling situation can, of course, also be treated in our 
Caldeira-Leggett like model. The main difference is that we need to expand our various 
fields not in the eigenstates of the trapping potential Xa( x ) but in the eigenstates x«( x ) °f 
the nonlocal Schrodinger equation 
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e' y' 



h 2 v 2 

2m 



+ V»(x)W(x) + J dx'Re[^( + )(x,x' ;e ^- / i)] X ' Q (x') (92) 



where e' a are the new eigenvalues and fi£( + )(x, x'; e) = J2 a ,a' Xa( x )^ia'( e )x^( x ')- In this 
new basis the nondiagonal elements of the selfenergies can now be neglected and we find 
essentially the same results as before. We only need to replace e a + S a by e' a and £ a (k) 
by J2 a ' (I <^ x Xa( x ')Xa'( x )) £a'(k). Neglecting the nondiagonal elements in this basis only 
requires that the real part of the retarded selfenergy is much larger that its imaginary part, 
which is always fulfilled in our case because k c A ^> 1. 

Summarizing, the coherent and incoherent dynamics of the gas in the trap is for our toy 
model quite generally solved by 

(<P«) (t) = (<f>a) (0)e-^-^V r ^ 2 (93) 

and 

N a (t) = N a (0)e~ rat + N a (l- e'^) , (94) 

respectively. In the limit £ — > oo, the average of the annihilation operators (<p a )(t) thus 
always vanishes but the average occupation numbers N a (t) relax to the equilibrium distri- 
bution N a = (ePte*-ri - l) -1 . Although this appears to be an immediately obvious result, 
its importance stems from the fact that it is also true if we tune the potential energy bias 
AV ex such that at low temperatures the groundstate Xg( x ) acquires a macroscopic occu- 
pation, i.e., N g ^> 1. The gas therefore never shows a spontaneous breaking of the U(l) 
symmetry, in agreement with the notion that we are essentially dealing with an ideal Bose 



gas in the grand-canonical ensemble [ 62] . The reason of the absense of the spontaneous sym- 
metry breaking can also be understood from our stationary solution of the Fokker-Planck 
equation in Eq. flSBD, which shows that the probability distribution for \<f> g \ is proportional 



to the Boltzmann factor e - ^^™^'^ 9 ' 2 in the degenerate regime of interest. Hence, the cor- 
responding Landau free energy F(\<f) g \) = {t' g — n)\4> g \ 2 never shows an instability towards 
the formation of a nonzero average of \<p g \ due to the fact that e' g — \i can never become less 
then zero. Once we introduce interactions between the atoms in the gas, this picture fully 
changes. 



III. SEMICLASSICAL THEORY OF BOSE-EINSTEIN CONDENSATION 

In the present and the folowing section we turn our attention to an interacting but homo- 
geneous Bose gas. However, in Sec. [V] we generalize our results and include also the effects 
of an external trapping potential. There are several reason for considering the homogeneous 
case first. One important reason is that for a homogeneous gas we are allowed to take 
the thermodynamic limit in which Bose-Einstein condensation becomes a true second-order 
phase transition. We are then in fact studying the dynamics of a spontaneous symme- 
try breaking under the most ideal circumstances, which may nevertheless have important 
implications for cosmology p3|-p6| and also lead to a better understanding of the recent ex- 
periments performed with liquid 4 He ||67| . Furthermore, the homogeneous case is in a sense 
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also the most difficult one because thermalizing collisions cannot lead to a macroscopic oc- 



cupation of the one-particle ground state |p8j|69|. Using the language of quantum optics, we 



also cannot compare Bose-Einstein condensation of a homogeneous gas with a single-mode 
laser and we are essentially dealing with the (in principle infinitely many) multi-mode sit- 



uation. Since realistic experiments with atomic rubidium and sodium gases p6| , |28| , |7D| are 
usually also in a multi-mode regime, it even turns out that an understanding of the homoge- 
neous gas is important for magnetically trapped gases as well. This is less so for experiments 



with atomic lithium |]2~7 |, because then the inhomogeneity plays a much more fundamental 
role. 

In the semiclassical theory of Bose-Einstein condensation in a homogeneous, interacting 
gas, the main emphasis is on the study of the coherent dynamics of the gas |f71 |. The rea- 



son for this is that if Bose-Einstein condensation is possible at all, it has to be caused by 
coherent processes since incoherent, i.e., kinetic, processes cannot lead to Bose-Einstein con- 
densation in the thermodynamic limit. The reason why Bose-Einstein condensation cannot 
be achieved in this way is easily understood. From the appropriate quantum Boltzmann 
equation (cf. Eq. ([81])) we obtain that the production rate for the number of particles in the 
zero-momentum state is given by 

dN(0;t) CV nT/ „ ^ 

_i^2= (l + jy(0;f , (95) 
at r el 

where r e i is the average time between two elastic collisions and C is a constant of 0(1). For 
initial conditions close to equilibrium but with iV(0; t ) = 0, we can in first instance neglect 
the evolution of the noncondensed part of the gas, i.e., the time dependence of r eh and we 
find that 

N(0; t) = {e c{t - to)l ^ 1 - l) . (96) 

Hence, a macroscopic occupation of the one-particle ground state is only achieved at times 
t > to + 0(r e i ln(iV)), which indeed diverge if iV — > oo. For the initial conditions of interest 
this condition can also be written as t > to + 0{r e i ln(L/A c )), with L the (linear) system size 
and A c = {2nTi 2 /mkBT c ) l l 2 the thermal de Broglie wavelength of the atoms at the critical 
temperature T c . 

Physically, we thus envisage the dynamics of the gas in response to evaporative cooling in 
the following way. In the nondegenerate regime the gas just thermalizes by elastic collisions 
in such a manner that the occupation numbers iV(k; t) are in a good approximation at all 
times equal to a truncated equilibrium distribution [|72| , [73| . Since the temperature of the 
gas is decreasing, this distribution becomes more sharply peaked at low momenta and we 
enter at some point the degenerate regime where nA 3 = 0(1) but no condensate has been 
formed yet because we know that this is impossible by means of elastic collisions. At this 
stage the most important question to be answered is: Is it nevertheless possible that for a 
noncondensed but degenerate distribution function N(k;t), the energy e'(0;t) of the zero- 
momentum state becomes less than the instantaneous chemical potential fi(t) and the gas 
becomes unstable towards Bose-Einstein condensation? If the answer to this question is 
no, Bose-Einstein condensation is impossible and the study of the dynamics of the gas is 
restricted to the normal phase and relatively well understood |]7^-[76|. However, if the answer 
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is yes, the dynamics of the gas after the point of instability is clearly more complicated and 
deserves further investigation. 

In view of the central importance of the above question for the topic of this paper, we 
now first of all show in Sec. [Ill A| that the gas indeed develops the required instability for 
Bose-Einstein condensation. We then turn to the dynamics of the unstable gas. In Sec. [111B 



the dynamics of the phase transition is first studied in the semiclassical approximation, in 
which fluctuations in the order parameter are neglected. Due to this neglect of fluctuations, 
the time at which the instability occurs in the gas is equivalent to the time at which Bose- 
Einstein condensation takes place. In Sec. [TV], however, we include the fluctuations into our 
theory and discuss in detail how this affects the semiclassical picture. 



A. Instability of the Bose gas 

From our Caldeira-Leggett toy model we remember that the (renormalized) energies e' a 
and the corresponding eigenstates Xa( x ) can be determined from a nonlocal Schrodinger 
equation once we know the retarded selfenergy fr£( + )(x, x'; t — t'). The same is true for a 
homogeneous atomic Bose gas, only with this exception that the nonzero selfenergy is now 
due to the interatomic interactions and not due to the presence of a reservoir. In a sense 
an interacting gas plays also the role of its own reservoir. Moreover, the homogeneity of 
the gas leads to an important simplification because translational invariance requires that 
the retarded selfenergy is only a function of the relative distance x — x'. Therefore, the 
Schrodinger equation in Eq. fl92|) is solved by Xk( x ) = Xk( x ) and 

e'(k) = e(k) + Re[ft£ (+) (k; e'(k) - fi)] , (97) 

where from now on e(k) = h 2 ~k 2 /2m again and 

£« (x - x'; t - f) = [ / E<+> (k; e y(M-x')-(*-OA) . ( 98 ) 

Having arrived at this result, it is important to briefly remind ourselves what it implies 
physically. Mathematically, we have assumed that the gas has at the initial time to the 
occupation numbers N(k) that correspond to a truncated equilibrium distribution with 
chemical potential \i. Subsequently, we have taken the limit of t — > — oo. The energy 
e'(k) is therefore the energy of the one-particle state Xk( x ) after all the transients have died 
out. However, as is shown explicitly in Ref. ||71|| , the memory time in the selfenergy is for 
the same reason as in our Caldeira-Leggett toy model only of 0{h/ksT) and hence much 
smaller than all the other relevant time scales in the problem. As a consequence we are for 
all practical purposes allowed to consider e'(k) as the energy of the state Xk( x ), given all 
the occupation numbers iV(k). The time dependence of e'(k;£) can then be seen as solely 
due to the time dependence of the occupation numbers N(k; t). 

Keeping the last remark in mind, we now proceed to solve Eq. (0). To do that we 
of course need an expression for the retarded selfenergy of a weakly interacting Bose gas, 
which follows once we know the full selfenergy ft£(k; t, t') defined on the Schwinger-Keldysh 
contour C°°. Unfortunately, even for a dilute system this quantity cannot be calculated 
exactly and some approximation is called for. The approximation that we will make here 
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is the so-called many-body T-matrix approximation that has recently been reviewed in the 
literature The main motivation for this approximation is that due to the smallness of 

the gasparameter (na 3 ) 1//2 it is very unlikely for three or more particles to be within the range 
of the interaction and we only need to take account of all possible two-body processes taking 
place in the gas. It goes beyond the much used Popov approximation |2^J7^] by including also 
quite accurately the effect of the surrounding medium on the effective interaction between 
the atoms in the gas [50[8I |, which turns out to be of crucial importance for our purposes. 



Given this effective interaction V(k, k', K; t, t') for the scattering of two atoms that at 
time t' have the momenta H(K/2 ± k') and at time t the momenta h(K/2± k), respectively, 
the exact selfenergy obeys the Hartree-Fock-like relation 

r dk' 

ft£(k; t, = 2% J —— V(k - k', k - k', k + k'; t, t')G(k'; t', t) , (99) 



where the Green's function equals again Eq. (p3|) but with e(k) replaced by e'(k) to make 
the theory self consistent. The retarded selfenergy therefore obeys 

ft£«(k;t-t') =2i ( (v {+ \k - k\k - k\k + k! ;t - t')G< {k! ■ t' -t) 

J (flu) 6 ^ 

+ V< (k - k', k - k', k + k'; t - t')G^ (k'; t' - t)) , (100) 

leading to 

r rlk' 

AE< + >(k; e'(k) - = 2 J — V<+'(k - k', k - k', k + k'; e'(k) + e'(k') - 2^)Af(k') 

if we also introduce the notation e~ for the limiting procedure e—iO. Note that the two terms 
in the right-hand side correspond roughly speaking to the contributions to the selfenergy 
from scattering processes out and into the state Xk(x), respectively. This will be more 
clear in Sec. [IV| when we discuss the imaginairy part of this expression that describes the 
(incoherent) dynamics of the occupation numbers. To see the instability of the gas, however, 
we only need to consider the real part. 

Hence, our next task is to determine the effective interaction V^k, k', K; t, t') by con- 
sidering all possible two-body scattering processes. Because we are dealing with bosons, 
the effective interaction is a sum of a direct and an exchange term and can be written 
as V(k, k', K; t, t') = (T(k, k', K; t, t') + T(-k, k', K; t, t' ))/2 in terms of the many-body 
T-matrix that obeys the Lippmann-Schwinger equation p2| 

T(k, k', K; t, if) = V(k - k')8(t, t') 

+T I dt " I V ( k ~ k")G(K/2 + k"; t, t")G{K/2 - k"; t, t")T(k", k' , K; t", t') , (102) 

Ti Jc°° J (Ztc)^ 

with V(k — k') the Fourier transform of the interatomic interaction potential. By iterating 
this equation, we immediately see that the many-body T-matrix indeed sums all possible 
collisions between two particles. Moreover, the Green's functions G(K/2±k"; t, t") describe 
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the propagation of an atom with momentum h(K./2 ± k") from time t" to time t in the 
gas. Therefore, we also see that the many-body T-matrix incorporates the effect of the 
surrounding gaseous medium on the propagation of the atoms between two collisions. This in 
contrast to the two-body T-matrix, that is well-known from scattering theory |83| and obeys 
the same Lippman-Schwinger equation but with the Green's functions G(K/2 ± k";t,t") 
replaced by 

G (K/2 ± k"; t, t") = -z e --( K / 2±k ")(*-*")/^ ( t) t ") 

describing the propagation of an atom in the absence of a surrounding gas, i.e., in a vacuum. 
Note also that the many-body T-matrix has a 5-function singularity on the Schwinger- 
Keldysh contour. As a result the retarded and advanced components are now given by the 
more general expression 

T (±) (k, k', K; t-t')= T 5 {k, k', K; t)5{t - t') 

±6(±(t - t')) (T>(k, k', K; t - t') - T<(k, k', K; t - t')) (103) 

but the Keldysh component still obeys 

T K {k, k', K; t-t') = (T > (k, k', K; t-t') + T<(k, k', K; t - t')) . (104) 

Before we turn to the solution of the Lippmann-Schwinger equation in Eq. (|102|) , it is 
instructive to briefly mention what we obtain if we now apply the pseudopotential method 



22fl . In this method it is assumed that the effective interaction equals 

V(k, k', K; t, t') = ^-5(t, t') , (105) 
m 

which implies that only the retarded and advanced components of the effective interaction 
are nonzero and both equal to 4nah 2 5(t — t')/m. Substituting this in Eq. ( |101| ), we simply 
find that 

/ \ „ . . 8imah 2 

^E (+) (k; e'(k) - /i) = (106) 

m 

and thus that e'(k;i) = e(k) + 8nn(t)ah 2 jm. On the basis of the pseudopotential method 
we therefore conclude that the Bose gas never shows an instability towards Bose-Einstein 
condensation because a momentum independent selfenergy only leads to a shift in the chem- 
ical potential, which plays no role in the dynamics of the gas, even if we are evaporatively 
cooling the gas and the density n(t) is a time- dependent quantity. For our purposes it is 
thus important to determine a more accurate expresion for the effective interaction. 

As a first step towards this goal, we consider the retarded and advanced parts of Eq. ( |102| ) . 
They are given by 

T( ± )(k, k', K; t-t') = Vik - k')5(t - t') 

+ ll dt " I ^ / ( k - k ")G' (±) (k ,/ ,K;t-t")T( ± )(k",k / ,K;t"-t') , (107) 
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where G^'(k",K;f — t") are the retarded and advanced components of the two-particle 
Green's function G(k", K; t, t") = G(K/2 + k"; t, t")G(K/2 - k"; t, t") describing the propa- 
gation of two (independent) atoms with momenta h(K/2 + k") and h(K/2 — k") from time 
t" to time t in the gas. Explicitly, we have that 

G^(k",K;t-t") 

= ±6(±(i - t")) (1 + iV(K/2 + k") + N(K/2 - k")) e -*(^0 £ ".K)-2M)(t-*")/» , (i 8) 

with e'(k", K) = e'(K/2 + k") +e'(K/2 — k") the appropriate two-particle energy. Performing 
a Fourier transformation we then obtain 

T (±) (k,k',K;e) = V(k - k') 

+ / ^ Wk - n l + ^(K/2 + k-) + iV(K/2-k-) (±) 

which is well-known in nuclear physics as the Bethe-Salpeter equation and was first used in 



the context of Bose-Einstein condensation by Brueckner and Sawada ||84|| . Compared with 
the usual Lippmann-Schwinger equation 

r rile" 1 

r(±>(k,k^) = y(k-kO + /g F v(k-k") e± _ L Hk/f) T^(k»,k>-,e) (no) 

for the two-body T matrix, we see that the Bethe-Salpeter equation incorporates the effect 
that if the intermediate states are already occupied it is more likely for bosons to scatter 
into these states. Indeed the factor 1 + N(K/2 + k") + iV(K/2 — k") arises as the net 
difference between two atoms scattering into and scattering out of the intermediate states 
XK/2+k"( x ) and XK/2-k"( x )> which are proportional to (1 + N(K/2 + k"))(l + N(K/2 - k")) 
and N(K/2 + k")A^(K/2 — k"), respectively. In a real space picture the effect of the Bose 
enhancement factors can be seen as a result of the fact that the scattering wavefunction of 
the two colliding atoms has to be matched to the many-body wavefunction of the gas at 
large interatomic distances fl85| . Seen in this manner, it will be clear that this effect is not 
small in the gasparameter (na 3 ) 1//2 and cannot be neglected beforehand if we aim to include 
all two-body processes in the gas. 

In general the solution of a Bethe-Salpeter equation is quite complicated. However, for 
the quantum gases of interest we can basically solve it analytically by making use of the 
fact in the quantum regime the temperatures are always such that the thermal de Broglie 
wavelength obeys a/A C 1. As a result the many-body T matrix at the relevant momenta 
and energies that are much smaller than h/a and h 2 /2ma 2 , respectively, can in a very good 
approximation be expressed as [77] 

rM(k,k<,K; £ ) = T 't ,(0 -° ; °> w (111) 

V ' ' ' 7 1 -T(±)(0,0;0)S(K;e) ' V ' 

where T^(0, 0; 0) = Anah 2 /m and we have also introduced the function 

E(K; f » = / Wf mK ' 2 + + N(K/2 - k " ] h-^J,K)-2») ' (U2) 
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Because H(K; e) contains only the principle value part of the integral, the above approx- 
imation gives only the real part of the many-body T matrix. This is all we need for the 
semiclassical theory, but for the full quantum theory we need also the imaginary part. The 
latter is found from the optical theorem | 83fl , which in the many-body case reads 

r dk" 

T«(k, k', K; e) - TH(k, k', K; e) = -2m j 5{e- (e'(k", K) - 2/x)) 

xT( + )(k,k",K; e )(l + iV(K/2 + k") + iV(K/2 - k"))T(-)(k", k', K; e) . (113) 

Note that no approximation is made to obtain this result. 

We are now in the position to evaluate the contribution of the V^G < -terra to the real 
part of the retarded selfenergy. Combining Eqs. ( |101| ) and ( |1 1 1| ) we obtain in first instance 

dk' T(+)(0,0;0) 



(2tt) 3 1 - T(+) (0, 0; 0)S(k + k'; e'(k) + e'(k') - 2/x) ' 

However, the dominant contribution to the integral comes from thermal momenta due to the 
presence of the distribution function jV(k'). It is not difficult to show that for these momenta 
rW(0,0;0)H(k + k'; e'(k) + e'(k') — 2/x) = 0(naA 2 ) and therefore, for temperatures near 
the critical temperature where nA 3 = 0(1), always much smaller than one. Under these 
conditions this contribution to the retarded selfenergy is essentially equal to 8Trnah 2 /m, i.e., 
the same as the full result obtained by the pseudopotential method. At this point we thus 
conclude that if the gas becomes unstable, it is due to the second term in the right-hand 
side of Eq. ( |101| ). This contribution we consider next. 

We first have to determine the component T < (k, k',K;e) of the many-body T matrix. 
In general Eq. ( p. 02] ) shows that 

T<>>(k, k', K; t - t') = % - J dt" J ^- V(k - k")G<>>(k" 7 K; t - t")T^\k" , k', K; t" - t') 
dt" J V(k-k")GM(k",K;t-t")T<'>(k",k',K;t"~t') , (114) 

which is immediately solved by 

T<'> (k, k', K; t - t') = % - j dt" J dt'" J (k, k", K; t - t" 



xG < ' > {k",K;t" -t'")T^\k",k',K;t'" -t') , (115) 



if we use the Bethe-Salpeter equation for the retarded component of the many-body T 
matrix. A Fourier transformation leads, therefore, in particular to 



r dk" 

T<(k, k', K; e) = -2ttz j — 5(e- (e'(k", K) - 2/i)) 



xrW(k, k", K; e)N(K/2 + k")N(K/2 - k")T ( -\k", k', K; e) (116) 

but also to 
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r dk" 

T*(k, k', K; e) = -27TZ J 8{e — (e'(k", K) — 2//))T^(k, k", K; e)T^(k", k', K; e) 

x (1 + iV(K/2 + k") + JV(K/2 - k") + 2N(K/2 + k")iV(K/2 - k")) , (117) 

which we need lateron in Sec. [TVj . 

The real part of the V^G^-term in the expression for the retarded selfenergy thus 
becomes 

V {+ \k, k', k") 2 N(k')N(k")- P 



r dk' 


r dk" 


1 (2tt) 3 J 


(2tt) 3 



e'(k) + e'(k' + k" - k) - e'(k') - e'(k") ' 

with V^(k, k', k") a convenient shorthand notation for the following matrix element of the 
effective interaction U (±) (k-(k'+k")/2, (k'-k")/2, k'+k"; e'(k / )+e / (k")-2/i). For momenta 
hk that are of the order of the thermal momenta h/A, this contribution to the retarded 
selfenergy is easily estimated to be only of O(2nT^ + \0, 0; 0)(naA 2 )) and therefore negligible 
compared tot the contribution from the V^G < -term that we considered previously and was 
shown to be equal to 2nT^(0, 0; 0). For thermal momenta we thus conclude that 

,,, . , 8im(t)ah 2 

e'(k; t) = e k + ^ . 118 

m 

However, this conclusion is not valid for momenta hk that are much smaller than the thermal 
momenta, because in that case the energy denominator in the integrant favors the small 
momenta where the occupation numbers are esspecially large in the degenerate regime. For 
such momenta, and in particular for the momenta obeying hk < h\/8irna <C h/A that are 
shown in Sec. [Ill B| to be the most important states for the dynamics of the gas after the 
occurrence of the instability, we find in a good approximation that 

,., , , 8irn(t)ah 2 

e'(k;t)=e(k) + 



m 



r dk' 


r dk" 


I (2tt) 3 J 


(2tt) 3 



yW(0, k-, k-)f N(k>; t)N(k"; t) ^ V . (119) 

Note that in the right-hand side V^ + \0,k',k") also implicily depends on t, because it de- 
pends on the distribution function iV(k; t) as shown for example in detail in Eq. QUI ). 



In principle, Eqs. ( |118|) and (|119|) already clearly show the tendency of the gas to become 
unstable towards Bose-Einstein condensation, because the energy of the one-particle ground 
state is shifted less upwards as compared to the one-particle states with thermal energies. 
To actually show when the gas is unstable we, however, need to compare the energy of the 
zero-momentum state with the instantaneous chemical potential. Using Eq. ( |11S| ) and our 
picture that during evaporative cooling the occupation numbers of the gas are almost equal 
to a truncated equilibrium distribution, the chemical potential becomes 

. , 8nn(t)ah 2 , , . 
M (t)~ ^ +fio(t) , (120) 

where the time dependence of the ideal gas chemical potential fio(t) = fio(n(t),T(t)) depends 
on the precise path in the density-temperature plane that is followed during the cooling 
process. An instability therefore occurs once the quantity 
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e'{0;t)-ii{t) ~ -fio(t)-2 



r dk' 


r dk" 


1 (2vr) 3 J 


(2tt) 3 



V«ro,k',k"} 



iV(k'; t)JV(k"; t) 



ft 2 k' 



k"/m 



71 that the 



becomes less than zero. It was the single most important conclusion of Ref. 
gas indeed develops the required instability for Bose-Einstein condensation if a > and the 
temperature is less than a critical temperature T c = Tq(1 + 0(a/Ao)) that is slightly higher 
than the critical temperature T of the ideal Bose gas. Although we can argue that this 
might be an artefact of the many-body T-matrix approximation, which certainly does not 
take proper account of all the critical fluctuations near the critical temperature, it should be 
noted that all the predictions of the many-body T-matrix approximation are fully confirmed 
by a renormalization group calculation that does take critical fluctuations into account and 
only neglects the effects of a mass renormalization |81fl . Moreover, the upward shift of the 
interacting critical temperature by an amount of O((a/A )T ) = (^((na 3 ) 1 / 3 ^) has recently 
also been found by path- integral Monte-Carlo simulations ||86|| . We therefore believe that the 
results obtained here give an, even more than qualitatively, correct picture of the behavior 
of the gas near the critical temperature. 



B. Semiclassical Dynamics 

Having positively answered the question if an interacting Bose gas can be experimentally 
quenched into the regime in which it is unstable for a spontaneous breaking of the U(l) 
symmetry, our next task is to consider the dynamics of the gas after it has developed 
the instability. At the semiclassical level our results obtained in Sec. [Ill Aj show that the 



effective action for the long- wavelength dynamics of the gas, i.e., for the states with momenta 
ftk < hy87ma, is given by 

S d [<j>\ <j>} = Jdt k (t) [ih^ - e'(k; t) + m*)) Mt) 

-^7/ E T (+) (O,O,O;O)^ /2+k (t)^ /2 _ k (t)0K/2-k'(t)0K/ 2+k '(t)l • (121) 

k,k',K J 

We can understand that this is the correct semiclassical action in two different ways. First, 
it is found from our Schwinger-Keldysh formalism if we realize that at the semiclassical level 
the 'path' followed by the gas on the backward branch of the Schwinger-Keldysh contour 
is identical to the 'path' followed on the forward branch. Applying this observation to 
the full effective action S e ^[ip* immediately leads to the action S cl [<f)*,<j)} describing the 
semiclassical 'path' on either branch. Second, we can make use of the fact that 0(x, t) can 
also be seen as the order parameter of the Bose gas. The corresponding time-dependent 
Landau-Ginzburg theory can easily be derived in the many-body T-matrix approximation, 
if we use the imaginary-time approach to the equilibrium theory of interacting Bose gases 
4j|52] . Peforming then the usual Wick rotation to real time, we again obtain S cl [(p*, 0] for 



the effective action of the order parameter. 

In particular, we thus find that the dynamics of the zero-momentum part of the order 
parameter, i.e., the condensate, is in first instance determined by 
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S C o l [^M =Jdt Ul{t) U^-e\Q-t)+v{t)\ Mt) ~ T(+)(O 2 ^° ;O) |0o(t)| 4 | • (122) 
Introducing the density no(i) of the condensate and its phase #o by means of the relation 



4>o(t) = J n (t)V e ld ° , this simply leads to 

SSKfi] = V Jdt (- e '(0;t)+ii(t) - r(+)( ° 2 0,0;0) ^o(t)) no® ■ ( 12 3) 

Notice that in the process of deriving the last equation, we have omitted the topological term 
/ dt (ihdno(t) / dt) that does not affect the equations of motion and is therefore irrelevant at 
the semiclassical level. It is in principle only important when we also want to consider the 
quantum fluctuations of the condensate. Note also that we have made use of the definition of 



the chemical potential, or equivalently of the Josephson relation [87], to take the condensate 



phase independent of time. If we do not measure our energies relative to the chemical 
potential, we must use <f>o(t) = yJn Q (t)Ve ieo ^ and UdOo/dt = —fJ>(t) instead, which then 
leads to exactly the same action. Clearly, the above action is minimized by no(t) = if 
e'(0; t) — fi(t) > 0. However, in agreement with our previous stability analysis, we have a 
nontrivial minimum at 

Mt) ~ "TW(0,0,0;0) ' (124) 

if e'(0; t) — n(t) < 0. This result gives the desired evolution of the condensate density after 
the gas has become unstable in terms of the time-dependent chemical potential //(£). Our 
next task is therefore to determine an equation of motion for the chemical potential. 

To achieve this we now also have to consider the interactions between the condensed and 
the noncondensed parts of the gas, which have been neglected thusfar. Substituting <pk(t) = 
0o(£)^k,o + </>k(£)(l — ^k,o) into the semiclassical action and integrating over the fluctuations 
0^ (t) describing the noncondensed part of the gas, we find that the exact semiclassical action 
for the condensate obeys 

S d [n Q ,n] = S d [n ,fi] - zh\n(Z d [n ,fi]) , (125) 

where Z cl [no,/j] represents the functional integral over the fluctuations for given evolutions 
of the condensate density and the chemical potential. Writing S c1 [4>q + 4>'*,4>o + 4>'} as 
Soi^ai 4>o\ + SfW* ■> 4>''i n o, lA we thus have explicitly that 



Z cl [n ,fi] = J <#'*]<#'] exp|-Sf[0'*,0';no,/4} • (126) 

With this action the total density of the gas is calculated in the thermodynamic limit as 

1 5S cl \n , u] r dk , . 

n < t )=^^r="°' i ) + /(W iV(k;t) ' (127) 

where the occupation numbers are found from 
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W(k; t) = ^ ^ = . (128) 



The latter two equations, together with Eq. ( |124| ), in principle give both the condensate den- 



sity rio(t) and the chemical potential (x(t) as a function of the total density n(t) and formally 
thus completely solve the semiclassical dynamics of the gas. However, the above derivation 
of the semiclassical theory appears to be not fully consistent. On the one hand we have in 
the calculation of the total density varied the full semiclassical action S cl [n , ji] with respect 
to the chemical potential, but on the other hand we have varied only S^lno, fj] with respect 
to no(t) to obtain the equation for the condensate. The following question therefore arises: 
Why does the term — ih\n(Z cl [no, fj]) not contribute to the dynamics of the condensate? 
Fortunately, the reason for this is well-known [S8|. Expanding — ihln(Z d [no, //]) in powers 



of no(t) we can easily convince ourselves that the fluctuations only lead to corrections to 
the linear and quadratic terms in Sq[tio, fj] that in fact have already been accounted for by 
using the renormalized energy e'(0;t) and the renormalized interaction T^ + ^(0, 0, 0; 0). To 
avoid a double counting of the effects of the interaction, we should therefore indeed neglect 
the fluctuations when calculating the condensate density. Formally, this is equivalent to the 
requirement that the condensate density fulfills the Hugenholtz-Pines theorem ]19| , because 



Eq. (|124|) can be rewritten as /i(t) = hH^\0;t) — hT^^^t) in terms of the usual normal 
and anomalous selfenergies in Sf[(p'*, (p';n , fi}. 

Determining the occupation numbers (0 k (i)0 k (£)) d requires solving an interacting quan- 
tum field theory, which cannot be done exactly. An approximation is thus called for. Taking 
only the quadratic terms in Sf[<p'*, <f)'; no, ft] into account amounts to the Bogoliubov ap- 
proximation. Indeed, the action for the fluctuations then becomes equal to 



S B [r,<P']= jdt jx>'k(*) ^-e(k)-n o (t)TW(O,O,O;O)^0L 



~T(+)(0, 0, 0; 0)no(t) £ (<j£(t)<j>\(t) + #_ k (t)&(t)) \ , (129) 



2 



if we use Eq. ( ]119p to evaluate the energy difference e'(k; t) — e'(0; t) = e(k) at the long 



wavelengths of interest here. The energies of the Bogoliubov quasiparticles thus obey 



huj(k; t) = \/e 2 {k) + 2n (t)T(+)(0, 0, 0; O)e(k) , (130) 

and are purely real at this level of approximation. However, in the Bogoliubov approxi- 
mation the quasiparticles are noninteracting. This is reasonable sufficiently far below the 
critical temperature when the condensate density is large, but not very close to the critical 
temperature |^1[ . In that case the interactions between the quasiparticles are very important 



and cannot be neglected. 

A qualitative understanding of the effect of these interactions can be obtained by making 
use of a simple mean-field theory. Since the condensate density is relatively small we can 
approximate the average interaction between the quasiparticles by 

1^7 E ^ (+ HO,O,O;O)0- /2+k (t)^ /2 _ k (t)0^ /2 _ k ,(t)^ /2+k ,(t)\ 

\ ZV k,k',K^0 / 
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~ 7rW(0,0,0;0) (n 2 (t) - n(t)n (t) - 7io(f)l £<#(0&(0H ■ ( 131 ) 

Within the context of a mean-field approximation, we therefore see that the interactions 
effectively renormalize the Bogoliubov theory to 

SrbW*^'] = Jdt jg 4£(t) {in^ -e(k)j ^(t) 

-irW(0, 0, 0; 0)no(t) £ (^W* k (*) + #-k(*)&(*)) 1 (132) 
and lead to quasiparticle energies 



nw(k; t) = Ve 2 (k) - (n (t)T(+)(0, 0, 0; 0)) 2 (133) 

that are purely imaginary at long wavelengths. The occupation numbers for these momen- 
tum states thus decay in favor of the condensate density. The typical time scale for the 
decay is of O(h/n T^ + \0, 0; 0)) and always much slower than h/k B T c because a/A c 1 for 
the quantum gases of interest. This a posteriori justifies our neglect of the memory effects 
in the above. In summary, the semiclassical picture of Bose-Einstein condensation can then 
be visualized as in Fig. |], where we have also tried to indicate the dynamics of the phase 
6o(t) due to the Josephson relation. After the condensate density has been formed in this 
fully coherent manner, the distribution function AT(k; t) is clearly not in equilibrium since 
all the particles in the low-energy states have been transferred to the condensate. In the 
third and last stage of the condensation process this distribution function relaxes to ther- 
mal equilibrium due to elastic collisions between the Bogoliubov quasiparticles. This can 
again be accurately described by a quantum Boltzmann equation |89|,p0| . We do not discuss 



this third stage in any detail here (see for example Refs. [[Fg,|K|]) and just end with a brief 



overview of the semiclassical picture of Bose-Einstein condensation. 



C. Overview 

In the semiclassical theory the divergence of the time required for Bose-Einstein conden- 
sation is avoided, because the incoherent dynamics of the occupation numbers N(k; t) for the 
one-particle states is coupled to the coherent dynamics of the energies e'(k; t) of these states, 
and the fluctuations of the condensate are neglected. As a result the process of Bose-Einstein 
condensation now proceeds in three stages. In the first stage the gas moves, for example as 
a response to quickly removing the hottest atoms, towards the critical region of the phase 
transition. This part of the evolution can be accurately described by a quantum Boltzmann 
equation, in which in a selfconsistent manner collisions between the atoms lead to changes 
in the occupation numbers AT(k; t) and these changes in their turn affect the mean- field 
corrections to the energies e'(k; t). As long as the gas is outside the critical region, however, 
these mean-field corrections are constant and we simply have e'(k;£) = e(k) + Sixnah 2 jm. 
Hence, if initially not enough energy is removed, the gas always remains in this regime and 
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then typically evolves towards equilibrium on a time scale of 0(r e i), which near the critical 
temperature is of 0((A c /a) 2 h/k B T c ) . 



If on the other hand sufficient energy is removed initially, the gas enters the critical 
region and the second stage towards Bose-Einstein condensation begins. In this second 
stage, the gas first develops the required instability for the phase transition. Physically, the 
instability is a result of the fact that in the critical region, fluctuations tend to decrease the 
mean-field corrections to the energies of the low-momentum states with fik < 0(h\/&Tma) 
80| , |81~|1 , which makes it possible for the zero-momentum state to acquire an energy which 



is less than the instantaneous chemical potential. It is then clearly energetically favorable 
to put a macroscopic number of particles in that state. It is important to note that for 
this instability mechanism, it is crucial that the typical time scale on which the mean-field 
corrections are established is only of 0(h/kBT c ) and therefore much faster than the time 
scale on which elastic collisions can change the occupations numbers. In that sense it can 



be said that Bose-Einstein condensation is 'nucleated' on a time scale of 0(h/kBT c ) [52 



As we have seen above, incoherent collisions alone cannot lead to a macroscopic occupa- 
tion of the zero-momentum state. After the gas has become unstable, the actual growth of 
the condensate density is therefore in the semiclassical theory caused by a coherent depletion 
of the low-lying excited states in favor of the one-particle ground state. At the end of the 
second stage, which also has a duration of 0((A c /a) 2 h/k B T c ), the gas has thus acquired a 
highly nonequilibrium energy distribution and has to come to equilibrium in a third and 
final stage. This last stage of the condensation process is again of a kinetic nature and can 
be accurately described by a quantum Boltzmann equation for the Bogoliubov quasiparticles 
89| , |90|1 . Due to the linear dispersion of these quasiparticles at long wavelengths the typical 



time scale for relaxation towards equilibrium is exceptionally long and of 0((A c /a) 3 h/kBT c 
as was shown by Eckern in the context of spin-polarized atomic hydrogen [^DJ . 

For completeness sake, we should mention that Kagan, Svistunov, and Shlyapnikov have 
put forward a different semiclassical picture and assert that during the second, coherent 
stage of the evolution, the gas does not form a condensate but, just as in a two dimensional 
gas, only a quasicondensate ||91|| . Mathematically, this implies that the gas does not acquire 
a nonzero average of the annihilation operator of the zero- momentum state ipo{t) alone, but 
only of the sum J2k<k ^k(£)e ik " x , where hk is a momentum cut-off of 0(h\f&ima) f24j] . This 
semiclassical picture is however not confirmed by our microscopic calculation, because from 
the above discussion of the instability mechanism it is clear that a possible instability for a 
quasicondensate is always preceeded by an instability for a real condensate since a removal 
of the long-wavelength fluctuations with k < ko always increases the mean-field corrections 
to the energies e'(k; t). 

The semiclassical theory neglects in particular the effect of thermal fluctuations on the 
condensate dynamics, which are nevertheless anticipated to be very important for conditions 
near the critical temperature when a large fraction of the gas is still uncondensed. For 
a quantitative comparison with experiments that are presently being performed [92], the 
inclusion of these thermal fluctuations is therefore clearly of interest. There are also several 
theoretical reasons why we need to consider the effect of thermal fluctuations. Although the 
semiclassical theory of Bose-Einstein condensation gives important insight into the dynamics 
of this phase transition, it has two important drawbacks. First, it makes a sharp distinction 
between the kinetic and coherent stages of the evolution, which in principle of course overlap 
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with each other. Second, the discussion of the growth of the condensate actually makes use 
of the presence of fluctuations, without taking these explicitly into account. This is due 
to the fact that at the semiclassical level we always have a solution to the equations of 
motion in which the gas remains at all times in the metastable, uncondensed phase. By 
considering the effect of thermal fluctuations, both these problems can indeed be resolved 
and we finally arrive at a single Fokker-Planck equation that describes all the three stages 
of the condensation process simultaneously. How this is achieved is reviewed next. 



IV. QUANTUM THEORY OF BOSE-EINSTEIN CONDENSATION 

In the semiclassical theory of the previous section the main emphasis is only on the 
average (4> )(t), whereas the aim of the quantum theory is to consider also the fluctuations 
and determine the full probability distribution P [</>*, 0; £]. Combining the results of Sees. |I| 
and |TJ, we find that this probability distribution can be expressed as a functional integral 
with an effective action that in the many-body T-matrix approximation is given by 



Mt") 



S e "ir,i>}= J ct d t'J ct dt" |E^ k (0 ^h^-e(k)+^t')y(t',t")-hE(k;t',t") 

~W E ^(k,k^K;^^^'>K/2 +k (O^K/ 2 -k(^0^K/ 2 -k'(t / 0^K/ 2+k '(^ // )l • (134) 

k,k',K J 

We next have to project the field ^(x, t) onto the real axis by means of t±) = </>(x, t) ± 
£(x, t)/2. Up to quadratic terms in the fluctuations £(x, t) and neglecting again the memory 
effects, which is also known as the Markovian approximation [f45| , [58| and is justified here 
because the memory time is at most of Oih/ksT), we then easily obtain 

S e ^[0*, 0; T, £] = £ f ^ dt' | - e ( k ) + MO - ftE<->(k)J 5 ky 

+h e ^ H (k, k', k")^ +k ,_ k (t')0k"(t')Uk'(t') 

v k" J 
+ E dt' { ~ e(k) + M0 ~ ^ S(+) ( k )) 



+h E ^ (+) (k, k', k> k , +k „_ k (t')<Mm MO 

v k" J 



^k,k' J -°° I 



+|E^(k,k',k")0 k , +k »_ k (t')0k"(OKk'(O , (135) 
^ k" J 

introducing the shorthand notation S < - ± - ) '^(k) for the quantities E^'^k; e'(k) — fi) that, 
just like the matrix elements V^ ,K (k, k', k"), implicitly also depend on time due to the 
time dependence of the occupation numbers. At this point it is possibly not immediately 
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clear why we do not need to consider cubic and quartic terms in the fluctuations. The reason 
for that has been explicitly mentioned in Ref. [[Fl[] and boils down to the general problem 
of working with an effective theory, i.e., we must be careful not to double count the effects 
of the interaction. This almost happened in Sec. [IllB] when we derived the semiclassical 
equation of motion for the condensate density and it would also happen here if we do not 
restrict ourselves to the linear and quadratic terms in the fluctuations. 

As long as all the one-particle states are stable, we can in first instance neglect the 
nonlinear terms in S e ^[(f)*, (p; £*, £] that are proportional to V^'^k, k', k"). Physically, 
this implies that we can treat the interacting atomic Bose gas as an (almost) noninteracting 
gas of dressed quasiparticles, 'Bose liquid'. As a result the effective action becomes 

essentially identical to that of our Caldeira-Leggett toy model in Sec. [II B| and we can 
immediately write down the corresponding Fokker-Planck equation. Denoting the imaginary 
part of ^S^(k) again by =F-R(k;t), we have 



ih^P[<f>*,<f>;t} 







£ — (e'(k; t) - iR(k; t) - /z(t))0 k P[0*, 0; t] 



d 



E 



d 2 



a£*(k)P[0*,0;t] . 



(136) 



However, to extract the physics from this Fokker-Planck equation, we now need to determine 
the quantities R(k;t) and HTi (k). Using our results from the semiclassical theory, in 
particular Eqs. (|101[) , ( \L13j) and ( p.l6| ), it is not difficult to show that the imaginairy part of 
the retarded and advanced selfenergies obeys 

r dk' r dk" 

R(k; t) = ^j^ysj^y^ t) + e'(k"; t) - e'(k' + k" - k; t) - e'(k; t)) 

x ^ (+) (k,k',k /, )| 2 [(l + iV(k , ;t) + Ar(k /, ;t))Ar(k' + k , '-k;t)-iV(k / ;t)Ar(k' / ;t)] , (137) 

which has recently also been obtained by an equation-of-motion method instead of the 
Schwinger-Keldysh formalism that is employed here [93|. To obtain the Keldysh component 
of the selfenergy we first have to realize that from the general Hartree-Fock-like relation for 
the selfenergy ft£(k;t, t') on the Schwinger-Keldysh contour it follows that 

' V K {k -k',k- k', k + k'; e'(k) + e'(k') - 2fi)N(k') 



hZ K (k; e'(k) 



+2 



//) = 2 

dk! 
(27p 



(2tt) 5 

V<ik - k', k - k', k + k'; e'(k) + e'(k') - 2//)(l + 2iV(k')) , 



138) 



after a Fourier transformation of the time difference t — t' . Subsituting now Eqs. ( 116 ) and 
( p.17 ) we then finally obtain 



/dk 1 c rlk" 
C2^ J C2^ 6 (e ' (k ' ; t] + e ' (k " ; t] ~ e ' (k ' + k " " k; t] ~ C ' (k; t)] 



dk" 



(2tt) 3 
1/W(k,k',k" 



(2vr 
[(1 



N(k'; t) + N(k"; t))N(k' + k" - k; t) 

+ Nik!- t)N{k"; t)(l + 2iV(k' + k" 



k;t))] . (139) 
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What does the above imply for the dynamics of the gas? Considering first the equation 
of motion of {(f>k)(t), the Fokker-Planck equation tells us that 

ih^ifaKt) = (e'(k;t)-iR(k;t)- f ,(t)){<j> k )(t) . (140) 

Hence the rate of decay of a particle with momentum Tik due to collisions with other particles 
in the gas equals T(k;t) = 2R(k;t)/h. This indeed precisely agrees with a Fermi's Golden 
Rule calculation, because the rate for a particle to scatter out of the momentum state Xk( x ) 
is 

4.7T r rile 1 r (71c" 

r-(k; t) = - I -— I -— 5 (e'(k'; t) + e'(k"; t) - e'(k' + k" - k; t) - e'(k; t)) 



x 



\/W(k,k',k") 2 (l + iV(k';t))(l + iV(k";t))iV(k' + k"-k;t) , (141) 



whereas the rate to scatter into that state is 

T"(k; t) = f I ^ / ^ 5 (e'(k'; t) + e'(k"; t) - e'(k> + k" - k; t) - e'(k; t)) 

x |\/ (+) (k,k , ,k")| 2 iV(k';t)iV(k";t)(l + Ar(k , + k ,, -k;t)) . (142) 

Note that both the in and the out rate are a factor of two larger than might have been 
expected in first instance due to the fact that the direct and exchange contributions interfere 
constructively in the collision process. 

Although the value of the decay rate T(k;t) = r ou *(k;t) — T m (k;t) is clearly a correct 
consequence of our Fokker-Planck equation, what is missing at this point is the dynamics of 
the occupation numbers N(k; t) themselves. This omission can be repaired by considering 
also the average (|0k| 2 )(^) = N(k;t) + 1/2. For that quantity the Fokker-Planck equation 
gives 

^(I0k| 2 >(t) = -2^(k;t)(|0 k | 2 >(t) - ^E^(k) . (143) 
As expected, this leads to the appropriate quantum Boltzmann equation 

|iV(k; t) = -(r^(k; t) - T m (k; t)) (jV(k; t) + \) + \iX K (k) 

= ~T mt {k; t)N{k; t) + T m (k; t)(l + A^(k; t)) , (144) 

if we make use of the observation that iL K {k) = T out (k;t) + r m (k;t). The gas, therefore, 
relaxes for t — > oo to an equilibrium in which the occupation numbers N(k) are equal 
to the Bose distribution evaluated at e'(k) — /x, where both e'(k) and \x are determined 
selfconsistently by the retarded selfenergy in the same way as in the semiclassical theory. 
An identical conclusion is reached by noting that in equilibrium the fluctuation-dissipation 
theorem is satisfied, i.e., ^S x (k) = — 2z(l + 2N(k))R(k). Therefore, the stationary solution 
of the Fokker-Planck equation equals 

^••*~l = n7v(i^ ex 4-7v(k)TT75W 2 }- < 145 > 
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In the degenerate regime and at long wavelengths the probability distribution for |0k| is 
thus proportional to the Boltzmann factor e _/3( - e ( k )-wlw _ Once again, this shows that the 
Landau free energy for the zero-momentum state F(\<p \) = (e'(0) — /i)|0 o | 2 is unstable 
when e'(0) < /i. If that happens we are no longer allowed to neglect the nonlinear terms in 
S e ff[(f)*, (ft; £*, £], because they are essential to ultimately stabilize the large fluctuations in 
the amplitude of 0(x). 

The linearized version of the full quantum theory given by S e ^[(f)*, <p] £*, £] thus ac- 
curately describes the first kinetic stage of the condensation process in which the gas is 
quenched into the unstable region of the phase diagram. The next question is: How does 
the quantum theory deal with the subsequent coherent stage? To answer that question we 
can neglect all the imaginary parts in the effective action, which are associated with the 
kinetic processes as we have seen, but now retain the nonlinear terms. The Fokker-Planck 
equation for the long-wavelength dynamics of interest then reads 

'(M -^))«/, k + I £ r(+)(o,o,o ; o)^ +k „_ k ^^ P[</> m ,M 
v k',k" / 

+ Ei (V(M-MWk+^ E rW(O,O,O;0)ft^^ k P[<P\M . (146) 

k 0( pk \ v k',k" / 

Since it contains only streaming terms the solution is immediately seen to be 

P[0*,0;t] = / P[|0i| 2 ;ti]II^(l^-0k'W| 2 ) , (147) 
J k 

where P[|0i| 2 ; ti] is the probability distribution after the kinetic stage, d (x, t) is the solution 
of the nonlinear Schrodinger equation 

in W^A = + 5(0; t) - /*(*) + T( + )(0, 0, 0; 0)|^(x, t)\ 2 ^j d (x, t) (148) 

with the initial condition d (x, ti) = 0i(x), and S(k; t) denotes the real part of the retarded 
selfenergy that is well approximated by the constant S(0; t) for small momenta Tik <C h/A 
and given explicitly in Eq. ( |119|) . 

The behavior of the solution in Eq. ( |147| ) has recently been investigated numerically by 
Damle, Majumdar, and Sachdev ||94|| . In agreement with the semiclassical picture, they 
find that the nonlinear Schrodinger equation has the property of depleting the states with 
momenta hk < h^87moa in favor of the growth of the condensate. In contrast with the 
semiclassical theory, however, the time it actually takes to form the condensate in a truly 
infinite system is infinite also. More precisely, the time to form a condensate in a system 
with linear system size L diverges as 

n ( lY 



n T(+)(0,0;0) \^ 
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introducing the correlation length of the condensate £ = 1/ 'y/8¥noa and the usual dynamical 
exponent z |)5| that was shown to be equal to 1 in this case [96]. From this outcome of the 



quantum theory we conclude that the fluctuations with wavelengths larger than the corre- 
lation length £ play a crucial role. In the semiclassical theory the effect of these fluctuations 
is only accounted for in a mean-field sense, which allows the condensation to take place in a 
finite time. Put differently, not taking fully account of the fluctuations on the longest length 
scales has the same effect as considering a finite-size system. For a truly infinite system the 
thermodynamic limit in fact prevents the formation of the condensate even though the gas 
becomes unstable in a finite amount of time. 

To understand the last remarks in a somewhat different manner, let us now consider the 
gas in a box with volume V = L 3 . Because of the energy gap of (2it?i) 2 / 2mL 2 between the 
one-particle ground state and the first excited state, there is always a regime (the details are 
given below) in which it is allowed to single out the zero-momentum state. In that regime 
the Fokker-Planck equation following from the zero-momentum part of S e ^[(j)*, 0; £*, £] is 

ih^P^l O ; t] = (s(0; t) - tR(0; t) - //(*) + ^T«(0, 0, 0; O)|0 O | 2 ) 0oP[0o, 0o; t] 

+ — (5(0; t) + iR(0; t) - p(t) + -T«(0, 0, 0; O)|0 O | 2 ) <P* P[<Pl 0o; t] 



1 d 2 



fiS K (O)P[0*,0o;t] • (149) 
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Anticipating now that the condensate dynamics is much slower than the dynamics of the 
noncondensed part of the gas, we can take for the selfeneries hT,^' K (0) and the chemical 
potential their equilibrium values. There is, however, again a slight complication because the 
streaming terms in the right-hand side of Eq. ( |149 ) show that the energy of the condensate is 
e'(0)+T(+)(0, 0, 0; O)|0o| 2 /^- To be fully consistent we must therefore evaluate frE^^O; e) 
at these energies to obtain the correct equilibrium values. Making then also use of the 
fluctuation-dissipation theorem that now becomes 

iR(0) = -4x^(0) (V(0) - fi + irW(0, 0, 0; O)|0 O | 2 ) (150) 

in the degenerate regime of interest, we find that the Fokker-Planck equation for a U(l) 
invariant probability distribution P[|0 o |;i] is 

z^P[|0o|;t] = -^(O)A^ (O )-^ + lT( + )(O,O,O;O)|0o| 2 ) </> P[\M,t] 
- f (/(0) -fi+ ^T (+) (0, 0, 0; O)|0 O | 2 ) <f P[\M,t] 

-^>^ P[IMt] - (151) 
The stationary solution therefore obeys 

f f r(+)(o,o,o ; o) IJL |2 \ ., |2 1 



P[\4>o\; oo] oc exp I -(3 I e'(0) - (j, H ' ' — |0 O 
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and shows explicitly that the nonlinear terms in the effective action S e ^[4>*, 0; £*, £] have 
indeed stabilized the Landau free energy for the condensate F(\4> \) = (e'(0) — fi)\(fio\ 2 + 
T( + )(0, 0, 0; O)|0 O | 4 /2V A in the case that e'(0) < \i. Note that this result also shows why 
the nonlinearities can be neglected in the stable regime when e'(0) > /i, because the most 
important feature of this stationary solution is that if we consider the probability distribution 
for the condensate density n = V\(po\ 2 , i-e., 

P[n ; oo] oc exp < -f3V I e'(0) - n H — n 1 n > , 

it is for a large volume V very sharply peaked around zero if e'(0) > fi. On the other 
hand, if e'(0) < \x it is very sharply peaked around the semiclassical solution no = (fi — 
e'(0))/T( + )(0, 0, 0; 0). This represents the phenomenon of spontaneous symmetry breaking 
in the context of Bose-Einstein condensation. 

To extract the time scale for the probability distribution to relax to equilibrium, it 



is convenient to realize that Eq. ( 151 ) is equivalent to the well-known Fokker-Planck 
equation for a single-mode laser |[45| . Indeed, introducing the dimensionless time vari- 
able r = t(iS^(0)/8)(2/?TW(0,0,0;0)/\/) 1/2 and in addition the dimensionless quantities 
/ = |0 O | 2 (2/3T(+)(O,O,O;O)/V) 1 / 2 and a = 2fi(\i - e'(0))(V/2/3T«(0, 0, 0; 0)) 1 / 2 , that in 
the laser theory correspond respectively to the 'laser intensity' and the so-called 'pump 
parameter', it acquires the standard form 

^PH; r] = jj {2(1 - a)I + 41^} P[I- r] (152) 

that has been studied numerically by Risken P7J . From the linear part of this equation we 



observe that initially the average 'intensity' behaves as 



(/>(r)^(e 2 --l) . (153) 



Hence, (J)(t) is of the order of its equilibrium value, i.e., of 0(a), when r is of C(ln(a)/a). 
The time it takes to form a condensate is therefore at least of 0(r e i ln(L/A c )) and again 
diverges for L — ► oo. At this point it is important to remember that the same time scale 
for the formation of a condensate was obtained in Sec. [Ill] from the quantum Boltzmann 
equation. This is not surprising because by considering only the dynamics of the zero- 
momentum state, a condensate can only be formed by incoherent collisions. Physically this 
is, however, only a good approximation as long as the mean- field interaction n T^+\0, 0, 0; 0) 
due to the condensate is small compared to the energy splitting of the one-particle states in 
the box or if 

/2tt\ 2 h 2 
M-e'(0)« — — . 



L J 2m 

If this condition is not fulfilled, which is generally the case for a large system size, we are 
no longer allowed to single out the zero-momentum state since all the states with momenta 
hk < h^8irn a are now coupled by the mean-field interaction. As we have seen in Eq. (|147|) , 
and in an approximate sense also in the semiclassical theory, the condensate can then also 
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be formed in a fully coherent manner. In the language of laser theory this means that Bose- 
Einstein condensation is then analogous to a multi-mode laser instead of a single-mode laser 
described by Eq. ( |152| ). 

The competition between the coherent and incoherent processes trying to form a conden- 
sate is of course contained in the Fokker-Planck equation associated with the full effective 
action S e ^[(j)*, 0; £*, £] given in Eq. (|135|) . By performing a Hubbard- Stratonovich trans- 
formation or by integrating over the fluctuations and quantizing the resulting theory, we 



ultimately find that [98] 
iK^P[<f>*,</>;t] = 

-T,4r ( ( e ( k ) + &£< + >(k) " MWk + ^ E ^ (+) ( k - k '> k")0 k , +k ,,_ k <^0 k < ] P[<p, 0; t] 
k % ^ V k , k „ 

+ E 9 



k 



( (e(k) + aEH(k) - (i(t))<fl + ^ E ^-)(k,k',k")0 k '0k"0k' + k"-k] 

\ k',k" / 

4 E Uz K (k)5 Kk , + I E V*(k, k, k")0 k , +k „_ k k ^ P[0*, 0; t] • (154) 

1 k,k' a <P\L a <P\L' \ V k" / 

We have already seen in detail how this equation describes the first and second stages of 
Bose-Einstein condensation. What remains to be shown is that it in fact also describes 
the third kinetic stage in which the condensed and noncondensed parts of the gas come 
to equilibrium. To do so, we make the ansatz that the probability distribution P[0*, <fi;t] 
separates as Po[0o> <fio', t\P\W*, <$'■> A- As w * n ^ e c l ear shortly, this essentially amounts to 
a Hartree-Fock approximation. Substituting this ansatz in Eq. ( 154 ) and integrating over 
the field 0'( x ) that i s again associated with the noncondensed part of the gas, we recover 
Eq. ( |149| ) for the probability distribution of the condensate 

ih^P [<f>* ,<f>o;t] = (W+)(o) -n(t) + ^tM(O,O,O;O)|0 o | 2 ) <f> P [<f>* ,(f>o;t) 
+ ^7 (fiE (_, (0) - m + -T<+>(0, 0, 0; 0)|^ o | 2 ) «P„[«. <A)I t] 

4aJk fiEA ' miCfe "' (155) 

if we are careful not to double count the effects of the noncondensed part of the gas that are 
already included in the selfenergies hH^' K (0). Integrating over O instead, we also arrive 
at a Fokker-Planck equation for the probability distribution Px[(f)'*,(f)';i\ that reads 

- E ^TF (( e ( k ) + ^ (+) ( k ) + 2V«(k, k, 0)n (t) - n(t))<H) P 1 [<j>'*, 0'; t] 
k^o 0( Pk 

+ ((e(k) + KZMQl) + 2V^(k, k, 0)n (0 - MWkj A[f, 0'; t] 

k^O 0( Pk 

"l £ o^^ )+V ^ k,0)fl0(i) '*^* 1 ' (156) 
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From the real part of the streaming terms we observe that in the presence of a condensate 
the renormalized energies of the states Xk( x ) with k 7^ are well approximated by 



e"(k; t) = e(k) + S(k; t) + 2n (t)T (+) (0, 0, 0; 0) , (157) 

as expected from a Hartree-Fock approximation. Furthermore, the quantum Boltzmann 
equation that follows from the imaginary terms in the right-hand side of Eq. (156) is again 
of the form 

^N(k; t) = -T™\k; t)N(k; t) + T m (k; t)(l + iV(k; t)) , (158) 

but T out and T in are now given by Eqs. ( pH|) and ( p2|) with iV(k' + k" - k; £) replaced 
by iV(k' + k" - k;t) + n (t)(27r) 3 5(k' + k" - k) and e'(k;t) by e"(k;t) to account for the 
macroscopic occupation of the zero-momentum state. Apart from the use of the renormalized 
one-particle energies in the expressions for the rates for particles to scattering out of and into 
the momentum state Xk(x), this kinetic equation is identical to the one studied by Semikoz 
and Tkachev |76[ . 



In equilibrium the occupation numbers are thus equal to a Bose distribution evaluated at 
e(k) + S(k) + 2n T^ + \0, 0, 0; 0) — ji. However, for a large system size we have seen that the 
chemical potential is given by the Hugenholtz-Pines relation y. = S(k) + n T (+) (0, 0, 0; 0). 
Substituting this we obtain that the Bose distribution is in effect evaluated at fiui{k) = e(k) + 
T( + )(0, 0, 0; 0)no, which precisely corresponds to the high-momentum (or Hartree-Fock) limit 
of the expected Bogoliubov dispersion Tiu(k) = (e 2 (k) + 2n T(+)(0, 0, 0; O)e(k)) 1 / 2 . For the 
quantum gases of interest the thermal momenta at the critical temperature are always much 
larger than hy/Wrmoa and the Hartree-Fock approximation presented above is very accurate. 
It can thus safely be used to describe the third stage of the condensation process. For 
applications of Eq. ( |154j ) at temperatures very far below the critical temperature, it is in 
principle also of interest to knows how we can recover the Bogoliubov approximation from our 



Fokker-Planck equation [ 9§| . Such extreme low temperatures have experimentally, however, 
only been obtained because of the inhomogeneity of the magnetically trapped atomic gases. 
Therefore, we do not discuss the Bogoliubov approximation here and return to this point in 
the next section when we have generalized the quantum theory to the inhomogeneous case. 
For now we only mention that the derivation of the Bogoliubov theory is not difficult in 
principle but unnecessary for a homogeneous atomic Bose gas, because our previous ways of 
looking at the Fokker-Planck equation in Eq. ( |154j) have already shown that it can accurately 
describe all the coherent and incoherent processes during Bose-Einstein condensation in that 
case. 



V. BOSE-EINSTEIN CONDENSATION IN AN EXTERNAL POTENTIAL 

Up to now we have only considered atomic Bose gases with effectively repulsive inter- 
actions, i.e., a positive scattering length a. The reason for this restriction is that a ho- 
mogeneous gas with effectively attractive interactions never develops an instability towards 
Bose-Einstein condensation as we have seen [ 100 ]. Instead it develops an instability towards 



a BCS-like transition that for bosons is known as an Evans- Rashid transition ||101|| . More 
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precisely, this state of affairs depends on the strength of the three-body interaction, which 
in principle can stabilize a condensate if it is sufficiently repulsive. Using the renormaliza- 
tion group approach of Ref . |^TJ , this can easily be studied quantitatively and the result is 
shown here in Fig. |2|. It is important to mention that the renormalization group calculation 
predicts the Evans-Rashid transition to take place in the mechanically stable region of the 
phase diagram, in contrast to the mean- field treatment ||102|| . Note also that in general 
the three-body interaction must be anomalously large to be able to lead to Bose-Einstein 
condensation and this is therefore not expected to occur in realistic atomic gases with a 
negative scattering length. 

Although this essentially rules out the achievement of Bose-Einstein condensation in a 
homogeneous gas with a < 0, this is no longer true if the gas is confined in an external 
trapping potential, as was first argued by Sackett and Hulet ||103|| and subsequently shown 
convincingly by Ruprecht, Holland, Burnett, and Edwards [ |104j| . In fact, the quantum theory 
for Bose-Einstein condensation in a trapped gas with attractive interactions is even more 
simple than for a gas with repulsive interactions, because the average mean-field interaction 
n\T^(0, 0; 0)| can never become much larger than the energy splitting of the trap due to 
the intrinsic instability of the gas to collapse to a dense phase. We therefore treat this 
case first and then turn to the discussion of Bose-Einstein condensation in a harmonically 
trapped gas with repulsive interactions, where the average mean-field interaction is usually 
much larger than the energy splitting in the trap and we cannot neglect the effect of the 
interactions on the one-particle eigenstates in the external potential. These conditions are 
in particular realized in the recent experiment by H.-J. Miesner et al. that has resulted in 



the first quantitative data on the dynamics of condensate formation |92 



A. Weak-Coupling Limit 

For a trapped atomic gas the appropriate one-particle states are no longer the momentum 
states Xk( x ) but instead the eigenstates x»( x ) °f t ne external potential V^x). As a result 
the Fokker-Planck equation for the interacting gas will now quite generally be given by 

a,a> 0( Pa [ p^i J 

+ E jL 1 0« - + + E vQ^fyh \ P[<p, 0; t] 

a,a> G< Pa { fifi* J 

E -ftTfttf \ nE *,*' + E V a,P;a',p4*^P' \ P[<f>* , 05 *] • (159) 
A a,a' °<PaO(p a , { ^ J 

For this equation to be useful, however, we now need to determine the various selfenergies 
and interactions. To do so, let us for simplicity first consider the normal state of the gas 
which is usually in the weak-coupling limit I^S^'I \e a — e a '| f° r the experiments with 
atomic alkali gases that are of interest here. In this regime we can neglect the nondiagonal 
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elements of the selfenergies and also the nonlinear terms in the right-hand side of Eq. (|159|) . 
The Fokker-Planck equation then reduces to 

ih^P[<f>*,frt) = 

a % a OtPa 

'E^W,^], (160) 



2 



« 



with hE^' K = KE^' K (e' a — fx) and e' a = e a + S a . This equation looks identical to the 
one we derived for our Caldeira-Leggett toy model but describes in fact quite different 
physics because now the selfenergies are not due to an interaction with a reservoir but due 
to the interactions between the atoms themselves. This is, however, completely hidden in 
the selfenergies as we show now explicitly. Performing the same manipulations as in the 
homogeneous case, we find in the many-body T-matrix approximation first of all that 



^ = (*ffiH) - = 2 £ V^(e' a + e',- 2fx)N, 



V<8;aAe) _ n , (161) 



2nh a ^^e--{e' a + e' -2fi) 



and similarly that 



KE£ = 2 £ V« M (e' a + e' p - 2fx)N p + 2 £ V<^(e' a + e' p - 2fx)(l + 2N P ) . (162) 

P (3 

Our next task is therefore the evaluation of the various many-body T-matrix elements 
V^p.'^'p, (e) = {T^0'p,(e) +T^0'^,(e))/2 that occur in these expressions, which is easily 
carried out with our experience of the homogeneous gas and by noting that the Lippmann- 
Schwinger equation in Eq. (|102[ ) is in fact an operator equation that only needs to be 
expressed in a new basis. In terms of the appropriate matrix elements of the interatomic 
interaction, i.e., 

Va,P;a>,P> = J dx J dx! ^(^(x') V(x ~ x');W (x) X/3 , (*') , (163) 

the retarded and advanced components of the many-body T matrix thus obey the Bethe- 



Salpeter equation [93 



T (±) / \ Tr T r l+Na'f+Np" (±) 

1 a,/3;a> ,/3> \ e ) ~ V a,J3;a',fi' + 2^ V a,/3;a» ,/3» ~T f~, —, _ 9j s 1 a" ,p»;a> ,(3> K e ) U b4 J 

a " ft" e \ e a" "T e i3" Z W 

and the optical theorem 



^a,/3;a',/3'( e ) ^aj;a',l3'{ e ) 



-27TZ £ 5(e - (e' Q „ + - 2#^(e)(l + N a „ + N p ii)T^„. al ^{e) . (165) 

a", 13" 
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Finally, we also need 



= ~ 2m E 5 ( e " «*" + 4" " 2//))^ j ;Q ,, /3 »(e)iV a »^,,T^^ ;Q , /3 ,(6) (166) 

a",j3" 

and the Keldysh component 

t^A") = ~ 2m E ^ - (4- + 4" - 2/*)) 

a",/?" 

x^SU/^Xl + N a „ + Np. + 2N a ,,Np,)T^,, a ,^e) . (167) 

Notice that the appearance of a 5 function in the last three equations in practice never 
leads to any problems, because in the experiments with trapped atomic gases the number 
of particles is always so large that the critical temperature of the gas is much larger than 
the energy splitting of the trap and the discrete sum over states is well approximated by an 
integral over a continuous spectrum. 

Making use of the above relations we again find that for the states with thermal energies 
the renormalized energy is given by the pseudopotential result 

. 8irah 2 c */\ / \ / \ 

e a W = e « + — ^— J « X Xa( X M X > %a(*J , (168) 

with n(x, t) the density profile of the gas cloud. For the states near the one-particle ground 
state there is, however, a correction given by 

r{+) 



2E E Ky.^A^NA^it)- 



(3 a" \P" 



and V^. all ^,,{t) = Va'!p;a",/3"( e 'a" + e '(3" ~ 2/x(t)). In the weak-coupling limit, the renormal- 
ization of the one-particle energies plays by definition not such an important role and we 
are much more interested in the imaginary part of the retarded and advanced selfenergies. 
They are given by 

R a (t) = 2nJ2i: 5(6> a „(t) + fr(t) - e' a (t) - e' p {t)) 

(3 a", 13" 

x|^ ( J! Q ,,^(0| 2 [(l + NAt) + NpMNpit) - N*,(t)Np,(t)] . (169) 
Furthermore, the Keldysh component of the selfenergies equals 

^ = -4*5: E WAt) + - - 

(3 a", 13" 

x|KSj!a",/9"(0| 2 [(l + N a -(t) + NpMNpit) + N a „(t)N p „(t)(l + 2N p (t))] . (170) 
As a result the quantum Boltzmann equation for the gas becomes 

|aU0 = -K ut (t)N a (t) + r£(f)(i + N a (t)) , (i7i) 

with the rate to scatter out of the state x»( x ) given by 
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Att 

K ut (t) = -EE S(e> a „(t) + - e' a (t) - e',{t)) 

11 13 a", /3" 

x|C^^(t)| 2 (l + N a „(t))(l + N^(t))N p (t) (172) 

and the rate to scatter into this state by 
Att 

C(*) = 1- E E <*(4»(*) + - e «W - 



^ /3 «",/3" 



x|^ Q „^(t)i 2 iV Q »(t)iV^(t)(l + iV>(t)) , (173) 



as expected from Fermi's Golden Rule. 

The quantum Boltzmann equation found above, in combination with the expressions 
for the shift in the energy levels, fully describes the dynamics of the gas in the normal 
state. However, close to the critical temperature when the occupation numbers of the one- 
particle groundstate Xc?( x ) start to become large, we are no longer allowed to neglect the 
nonlinear terms in the Fokker-Planck equation in Eq. ( 159 ). Applying again the Hartree- 
Fock approximation, which is essentially exact in the weak-coupling limit as we have seen in 
the homogeneous case, and substituting the ansatz P[<f>*, (ft; t] = Po[<f> g , 4> g ', t}P\[(p'*, (ft'; t] into 
Eq. ( |159| ), we obtain for the probability distribution of the condensate 

ih^PoWg, 0,5 1] = ~ (e g + ffiW - A*(t) + T^t)^ 2 ) cj> g P [^ ^ t] 
+ A (e g + ksH _ ^t) + T^jno, 2 ) <j> g ; t] 

' ^-r-REfPoW,^;*], (174) 



2. 



where it should be remembered that for consistency reasons the selfenergies need to be 
evaluated at the energies e' (t) + T^> g (t)\(f)g\ 2 . In equilibrium the fluctuation-dissipation 
theorem then reads 

iR* = "I^f (4 - A* + V mi o„ 2 ) , (175) 

which guarantees that the probability distribution for the condensate relaxes to the correct 
equilibrium distribution 

P o [0*, 9 ; oo] oc exp {-/? (4 - n + 2|^|0 3 | 2 ) \<f>A . 

The Fokker-Planck equation, and hence the Boltzmann equation, for the noncondensed part 
of the gas is the same as in the normal state because for a discrete system the contributions 
from the condenstate to the selfenergies is already included in our expressions for hH^' K 
and needs not be separated out explicitly. This implies of course, that if we evaluate the 
rates r° ut (i) and T™(t) by replacing the sum over states by an integral over a continuum, 
we first need to separate the contribution from the condensate to obtain the correct result. 
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In summary we have thus found that in the weak-coupling limit the dynamics of Bose- 
Einstein condensation in an external trapping potential is described by the coupled equations 
in Eq. (|171 ) and Eq. ( |174 ), i.e., by a Boltzmann equation for the noncondensed part of the gas 



and a Fokker-Planck equation for the condensate, respectively. If we apply these equations 
to a Bose gas with a positive scattering length, we essentially arrive at the theory recently 
put forward by Gardiner and coworkers ||105|| . Indeed, neglecting the energy shifts and 
assuming that the noncondensed cloud is in equilibrium with a chemical potential fi > e g , 
our Fokker-Planck equation for a U(l) invariant probability distribution of the condensate 
becomes 

in^p.UM = -£ nsf A (e g - ai + t^ 9 ,M 2 ) <P 9 PoW;t} 

~ f^^S ( c * " a* + ri^l^l 2 ) 0^Po[|0J; t] - i^f^^^n^l; *1 • ( 176 ) 

It is again equivalent to the Fokker-Planck equation for a single-mode laser and thus incor- 
porates in a slightly different language exactly the same physics as the theory of Gardiner 
et al. A comparison with the experiments of H.-J. Miesner et al. shows that the theory of 
Gardiner and coworkers is unfortunately not fully satisfactory ||92|| . The reason for this is 
presumably that these experiments are really in the strong-coupling limit, where the above 
theory does not apply. As mentioned previously, the theory is however applicable to gases 
with a negative scattering length. Since we know from our experience with homogeneous 
gases that the chemical potential will never be larger than e' g in this case, we clearly see 
from the equilibrium probability distribution Po[<j)*,<j)g; oo] that a condensate in a gas with 
effectively attractive interactions is only metastable and will ultimately collapse to a dense 
state after it has overcome an energy barrier, either by thermal fluctuations or by quantum 
mechanical tunneling. In the next section we study the dynamics of this collapse by an 
extension of the method that has been devised in Ref. |106| and that is more accurate than 
the weak-coupling theory discussed sofar. A brief review of this approach also turns out to 
be a useful intermediate step towards the discussion of the strong-coupling limit presented 
in Sec. [VJj. 



B. Negative Scattering Length 

For defmiteness we now consider an external potential in the form of an isotropic har- 
monic oscillator. Thus V ex (x) = mo; 2 x 2 /2 and the weak-coupling limit then requires that 
8nh 2 \a\n(0,t) /m <C fvuj. During the growth of the condensate this condition is essentially 
always fulfilled but during the collapse, that occurs after the number of condensate particles 
has become too large, the density in the center of the trap quickly increases and we leave 
the weak-coupling regime. To study also the dynamics of the collapse requires therefore a 
more accurate theory. However, an important simplification occurs because the typical time 
scale for the collapse turns out to be of 0(1/ u) which is much shorter than the time between 
collisions r e \_. Indeed, their ratio is of order 

n(0,t)T(+)(0, 0; 0) / a \ 
hto \AJ 
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and the dynamics of the condensate is thus collisionless. From our general Fokker-Planck 
equation in Eq. (|159|) we then see that the condensate wavefunction obeys 



E 



aOW + M^i,}) (^>(t) + £ v^A^W^'W){^){t) ) ■ (177) 



In coordinate space this corresponds, in a reasonable first approximation, to the expected 
nonlinear Schrodinger equation at nonzero temperatures 



h 2 V 2 
2m 



rn 



n'(x,t) 



4nah 



rn 



|<0(x))(t)| 2 <0(x))(t) , (178) 



that was first studied by Goldman, Silvera, and Leggett ||107|| and Huse and Siggia ||108 
for a gas with positive scattering length, because these authors were at that time interested 
in spin-polarized atomic hydrogen. It was recently applied to atomic 7 Li by Houbiers and 
Stoof [|109|| and Bergeman ||110|| . Notice also that we have evaluated the sum over states by 
an integral over a continuum to obtain the noncondensate density n'(x, t) in the right-hand 
side of this equation. 



Neglecting the variation of the noncondensate density on the size of the condensate [109 



we conclude that the dynamics of the collapse is, apart from an unimportant shift in the 
chemical potential, determined by the Gross-Pitaevskii equation ||111 | 



^|(0(x))(t) 



h 2 V 2 A-Kah 2 1 

+ V"(x) -fi+ |(0(x))(t)| 2 (0(x)>(t) 



2m 



in 



(179) 



More precisely, this determines only the semiclassical dynamics. If we also want to study 
the quantum fluctuations, which is necessary if we are also interested in how the conden- 
sate tunnels through the macroscopic energy barrier, it is most convenient to calculate the 
grand canonical partition function of the condensate | 112 |. Quantizing the Gross-Pitaevskii 
equation we obtain for this partition function the functional integral 



Z gM = J d[ip*]d[i)] expj--^^*,^] 
over the complex field ^(x, r) and with the Euclidian action 

J dr J <ix ^*(x, r) 



180) 



9 h 2 V 2 Trex , , 2rf 

jr + ^7— + V ex x -fi+ 

or 2m m 



(x,r)| 2 U(x,r) . (181) 



As always for Bose systems, the integration is only over fields that are periodic on the 
imaginary time axis. 

Although it has recently been shown by Freire and Arovas that the tunneling process 
can also be studied in terms of the complex field ip(x., r) ||113| , we believe that it leads 
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to somewhat more physical insight if we use instead the fields p(x, r) and 0(x, r) that 
correspond to the density and phase fluctuations of the condensate, respectively. They are 
introduced by performing the canonical variable transformation [^4] 



^(x,r) = ^/p(x,r)e ie ^ 
in the functional integral for the partition function. As a result we find 

Z gr (fi) = J d[p}d[9] exp{-^[p,£;/i]} , (182) 



with 



S E [p,e;n} = /"drfdx ^p(x, r)^|^ + ^^(V^(x, r)) 2 

+ , r , (Vp(x,r)) 2 + ^(x)p(x,r)-pp(x i r) + ^p 2 (x,r)') . (183) 
8mp(x, T) m J 

Next, we notice that this action is only quadratic in the phase fluctuations. The field (9(x, r) 
can therefore be integrated over exactly, because it only involves the evaluation of a gaussian 
integral. 

Compared to ordinary gaussian integrals there is, however, one slight complication which 
is associated with the fact that #(x, r) are phase variables. This implies that the periodicity 
of the original field ^(x, r) only constraints the phase field #(x, r) to be periodic up to a 
multiple of 2n. To evaluate the grand canonical partition function in Eq. (|182[) we must 
therefore first integrate over all fields 0(x, r) that obey the boundary condition 0(x,h/3) = 
6*(x, 0) + 2irj and subsequently sum over all possible integers j. Because these different 
boundary conditions only affect the zero- momentum part of 9{x, r) we first have to evaluate 
the sum 

E J d[9 ] exp|-z ^ dr N (r)-^j , 

with Nq(t) = J dx p(x, r) the number of condensate particles. After performing a partial 
integration on the integral in the exponent, we can carry out the path integration over 9q(t) 
to obtain 



^2 e 2 ™N j$ 



dN (r) 
Or 



As expected, the integration over the global phase of the condensate leads to the constraint of 
a constant number of condensate particles, i.e., Nq(t) = Nq. Moreover, we have J2j e 2mN °^ = 
X),-5(iVo — j), which restricts the number of condensate particles to an integer. Putting all 
these results together, we see that the integration over the zero-momentum part of p(x, r) 
is only a sum over the number of condensate particles and we have that 

ZM=Y,e^ N °Z No . (184) 

N 
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Here we introduced the canonical partition function of the condensate, which is apparently 
equal to the functional integral 



Zn = J d[p]d{6] eX p{-^[p,0;O]} (185) 

over all the nonzero momentum components of the density and phase fields. 

The integration over the nonzero momentum components of the phase field #(x, r) is 
easily performed, because it now involves an ordinary gaussian integral. Introducing the 
Green's function for the phase fluctuations G(x, x'; p) by 

^ ((Vp) ■ V + P V 2 ) G(x, x'; p) = *(x - x') , (186) 
we immediately obtain the desired effective action for the density field 

8M = r*!* 1 !'" R^^W;,)^*' ^ 



lo J J \ 2 dr dr J 

+ 1 (^yW X ,r)f + nx W x,,) + — Ax,r)} . (187) 

Being an action for the density fluctuations of the condensate, Se[p] also describes all the 
collisionless modes of the condensate. This is important for our purposes, because the mode 
which becomes unstable first, determines precisely how the condensate collapses. Moreover, 
it determines the probability with which the collapse is going to take place, both for quantum 
and thermal fluctuations, since the energy barrier is smallest in that direction of the config- 
uration space. It should be noted that as long as we can neglect the interaction between the 
condensate and the thermal cloud, the action Se[p] describes also the collective modes of a 
gas with positive scattering length. For various other theoretical approaches that have been 
applied under these conditions see, for example, Refs. | 114 - 122|| . The actual measurements 



of the collective mode frequencies have been performed by Jin et al. ||123|| and Mewes et al. 



|1 2 M and are at sufficiently low temperatures indeed in good agreement with the theoretical 



predictions ||125| , |126|1 . We expect the same to be true for a gas with effectively attractive 



interactions and, therefore, the action Se[p] to be a good starting point for the following 
discussion. 

To obtain the collisionless modes explicitly we consider first the case of an ideal Bose gas 
by putting a = 0. For the ideal Bose gas we expect the gaussian profile 




p(x; q(r)) = N —— exp -— - (188) 



to describe an exact mode of the condensate. The reason is that in the noninteracting case we 
can make a density fluctuation by taking one particle from the condensate and putting that 
in one of the excited states of the external potential. The corresponding density fluctuation 
obeys 

-i{e a -e g )t/h^ 



5p(x,t)cxe-^-^ x ;(x) X9 (x) . 
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For our isotropic harmonic oscillator a combines the two angular momentum quantum num- 
bers I and m with the quantum number n that counts the number of nodes in the radial 
wavefunction Xnt{ x )- The density fluctuation then becomes 



5p(x,t)oce-( 2 " + ^^(x)r; m (x 



e 



-x 2 /2l 2 



(ttZ 2 ) 3 / 4 



with e n £ m — e oo = (2n + £)hto the excitation energy and I = (h/mto) 1 ^ 2 the size of the 
condensate wavefunction. Comparing this now with the expansion of the gaussian profile in 
Eq. (|188| ) around the groundstate density profile, which is obtained by substituting q{r) = 
I + 5q(r), we find that 



5a(r) e~ a "I 21 
<fp(x, r) = -V6N ^ Xl0 (x)Y* (x)-— 4 (189) 

has precisely the same form as a density fluctuation in which one particle is taken from the 
condensate and put into the harmonic oscillator state with quantum numbers {nlm) = (100). 
The frequency of the so-called 'breathing' mode described by the gaussian density profile 
must therefore be equal to 2u. 

To proof that this is indeed correct, we need to evaluate the effective action Se[p\, and 
hence the Green's function G(x, x';p), for a gaussian density profile. Substituting such a 
profile in Eq. (|186|) immediately leads to G(x, x'; p) = G(x, x'; g)/p(x'; q), with 

- (-- |x • V + V 2 ) G(x, x'; q) = 5(x - x') . (190) 
m \ q z J 

The latter equation can be solved if we can solve the eigenvalue problem 

( v2 - ?3 ^ ^« x > ■ < 191 > 

This turns out to be an easy task, because substituting 

e x 2 /2q 2 

£ n£m (x) = £ n£ (x) Y lm (±) (192) 

X 

gives essentially the radial Schrodinger equation for an isotropic harmonic oscillator with 
frequency uj q = h/mq 2 , i.e., 

2m ( h 2 d 2 1 2 2 h 2 m + l) 3 + \ _ , . . , , . 

" IF {- 2^dx~ 2 + 2 m ^ X + ~ 2^ ) U{X) = XneU{x) ■ (193) 

The desired eigenfunctions are therefore £ n «m( x > <l) — ( Pn^m(x)e x2/ ' 2 ' j2 , with <p„£ m (x) the usual 
(properly normalized) harmonic oscillator states with the energies (2n + i + 3/2)huj q , and 
the corresponding eigenvalues are X n e(q) = —2{2n + £)/q 2 . Introducing finally the 'dual' 
eigenfunctions £n&n( x ;<7) = l P*n£m{*) e ~ x ^ 9 > the Green's function G(x, x';g) is given by 

G(x, x'; q) = 6^m( x ! glrrrrUfx'; q) ■ (194) 
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Note that prime on the summation sign indicates that the sum is over all quantum numbers 
except [nlrn) = (000). The latter is excluded because the associated eigenfunction £000 ( x ; q) 
is just a constant and thus does not contribute to G(x, x'; p), which is defined as the Green's 
function for all phase fluctuations with nonvanishing momenta. 

Putting all these results together, we see that the dynamics of the collective variable q(r) 
is determined by the action 



S E [q] 



dr 



{3mN Uq 



lf / 3r 3 22 



Tu3 



dr 




+ V(q) 



(195) 



that is equivalent to the action of a particle with effective mass m* = 3mN Q /2 in a potential 
V(q) = N (3h 2 1 Amq 2 + 3muj 2 q 2 / A) . As expected from our previous remarks, this potential 
has a minimum for q = I and can be expanded near its minimum as 



V(q)-^N hu J +^m*(2u J ) 2 (6q) 2 



(196) 



It thus comfirms that the gaussian profile describes a breathing mode with frequency 2uj 
around an equilibrium density profile that is given by p(x; I) = A^ |xooo( x )| 2 . 

Up to now we have only considered the breathing mode of the condensate. However, 
the explicit expression for the Green's function G(x, x';g) in Eq. ( |194| ), together with our 
argument for the form of the density fluctuations in the various modes, presents us with a 
clue on how to obtain all the collisionless modes. It suggests that we should expand the 
density profile of the condensate as 



I \ I ( W 1 n I \£riim(x; <?(r)) 

p(x,r) =p(x;g(r)) + ^ C nlm {r) (7n?2(r)) 3/4 



(197) 



where the sum is now over all quantum numbers except (n£m) = (000) and iniva) = (100). 
Substituting this in the effective action Se[p] and expanding up to second order in the 
complex coefficients C n £ m , we find after somewhat tedious but straightforward algebra ||127 
that 



S E [p]^S E [q] + T," r dT JN 



mq 
2n + 



dC, 



nim 



dr 



+ 



mq A 



:{2n + £)\C n£m \ 



(198) 



Neglecting then the interaction between the breathing mode and all the other modes of 
the condensate, we conclude that the coefficients C n f m are harmonic oscillators with the 
frequencies 



u n£m = (2n + 1) 



h 

ml 7 



(2n + i)u . 



(199) 



Hence, we have indeed succeeded in correctly describing also all the other collective modes 
of the condensate. At this point, we might wonder why we have chosen in Eq. ( |197p to 
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still treat the breathing mode with the collective variable q(r) and not with an additional 
coefficient C wo (r) and the expansion 

p(x,r) = p(x;Z) + £'a^(T)^M (200) 

•aim \ I 

Doing that is clearly also a valid procedure and results at the quadratic level in 



n£m 



mZ 2 



4iV \ 2n + 



+ ^(2n + £)|C n ^| 2 ) • (201) 



The reason for using the former, more general, choice is related to the fact that when we 
include the attractive interaction, the breathing mode turns out to be the unstable mode of 
the condensate. It therefore plays a special role in the dynamics of the collapse, as we will 
see now. 

Considering first again only gaussian density profiles, the action Se[q] is again that of a 
particle with effective mass m* = 3miVo/2 but now in the potential [ |128| ] 



,i 3h 2 3 9 9 iVn h 2 \a\ , 

V(q)=N [- + -mw 2 q 2 - ° ' 1 . 202 

i 4mg 2 4 ^/2it mq 3 1 



The physically most important feature of this potential is that it is unbounded from below, 
since V(q) — > — oo if q j 0. Hence, the condensate always has the tendency to collapse to the 
high-density state lim^o p( x ; q) = N 5(x.). However, if the number of condensate particles 
is sufficiently small, or more precisely if [|129|[ 



*b < ^ A * 0.68 ^ , (203) 
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the condensate has to overcome a macroscopic energy barrier before it can collapse. Un- 
der these conditions the condensate is therefore really metastable and can in principle be 
observed experimentally. The most important question in this respect is of course: How 
metastable is the condensate? Within the gaussian approximation this question is easily 
answered, because then the dynamics of the condensate is equivalent to the dynamics of a 
particle in an unstable potential, as we have just seen. We therefore only need to evaluate 
the WKB-expression for the tunnneling rate | |130| | and compare this to the rate of decay 



due to thermal fluctuations by calculating also the height of the energy barrier. The out- 
come of this comparison for the conditions of the experiment with atomic 7 Li is presented in 



Ref. [131 1 and shows that for the relatively high temperatures T 3> Tiuj/kB that have been 
obtained thusfar |27j the decay by means of thermal fluctuations over the energy barrier is 
the dominant decay mechanism of the condensate. 

More important for our purposes, however, is that sufficiently close to the maximum 
number of condensate particles N m the collective decay of the condensate discussed above 
is always much more probable than the decay due to two and three-body collisions that 
lead to a spin-flip or the formation of 7 Li molecules, respectively. As a result the collapse 
of the condensate should be observable within the finite lifetime of the gas. In fact, on the 
basis of this separation of time scales we expect the condensate to go through a number of 
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growth and collapse cycles ||131| , |132j| . Theoretically this physical picture arises as follows. 
Starting from a gas with a number of atoms N 3> N m , the initial growth of the condensate 
as a response to evaporative cooling is described by the Boltzmann equation in Eq. (|171 ) 



and the Fokker-Planck equation for the condensate in Eq. (|174|) . However, if the number 



of condensate atoms starts to come close to N m , the condensate fluctuates over the energy 
barrier and collapses in a very short time of 0(1 /u) \\133 \. During the collapse the condensate 



density increases rapidly and two and three-body inelastic processes quickly remove almost 
all the atoms from the condensate. After this has occurred the condensate grows again from 
the noncondensed part of the gas and a new growth and collapse cycle begins. It is only after 
many of these cycles that enough atoms are removed for the gas to relax to an equilibrium 
with a number of condensate particles that is less than N m . This is shows quantitatively in 
Fig. [| for the experimental conditions of interest. 

A final issue which needs to be addressed at this point is the actual dynamics of the 
collapse and in particular how we must include the effect of the inelastic decay processes on 
this dynamics. From our Caldeira-Leggett toy model we know that the inclusion of inelas- 
tic processes in first instance means that we must replace the collisionless Gross-Pitaevskii 
equation for the condensate by a dissipative nonlinear Schrodinger equation with noise. Con- 
sidering for simplicity only two-body relaxation since three-body recombination processes 
can be treated analogously, we now have to deal with the stochastic equation 

{ m ¥t + yex(x) + " ~ (^T " lhG °) ^ (x ' t)|2 } = ^ x ' t} ' (204) 

where the strength of the noise r?(x, t) is such that 
dN (t) _ d r d 



dt 



= -Jdxn Q (x,t) = jjd* (|0(x)| 2 )(t) = -2G J dxng(x,t) . (205) 



The latter equation shows that the condensate density no(x, t) obeys the expected rate 
equation 

J^no(x,t) = -2G nl(x,t) (206) 

irrespective of the fact that for atomic 7 Li we have that Go ATr\a\h/m and hence that the 
equation for (</>(x))(i) is in a good approximation still equal to the Gross-Pitaevskii equation 
in Eq. ( |179|) . On the basis of these arguments, which are in fact quite similar to the ones 
used for the strong-coupling limit of our Caldeira-Leggett toy model, we then conclude that 
the dynamics of the collapse can still be studied by the effective action Se[p\- We only need 
to add the constraint that 



dN (t) 



dt 

In the gaussian approximation this simply becomes 



2G /dxp 2 (x,t) . (207) 
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which needs to be solved in combination with the dynamics of the condensate 'width' q(t) 
obeying Newton's law 

*d 2 q(t) dV , , ., . 

m "ST = ' (2M) 

A numerical solution of these two coupled equations for N (t) and q(t) shows that, in con- 
trast to the case when there is no decay ||133|| , the collapse is not complete and is finally 



arrested when the number of condensate particles becomes too small ||131|| . In the gaussian 
approximation we find that the final number of condensate atoms is of 0(1). Preliminary 
experimental results indicate that this number is presumably too small and about 1% of the 



condensate atoms remain after a single collapse [134]. This difference between theory and 



experiment is presumably due to the fact that the gaussian approximation to the density 
profile of the condensate overestimates the densities in the final stage of the collapse when 



the decay of the condensate is most severe [135 



Fortunately, our discussion of the ideal Bose gas shows how one can systematically im- 
prove on the gaussian approximation by including also the effect of the other collective 
variables C n £ m (t) in the density expansion in Eq. ( |197| ). Of course, if one is only interested 
in the dynamics of the collapse it is also possible to (numerically) solve the Gross-Pitaevskii 
equation together with 



dN (t) 
dt 



-2G |rfx|(0(x))(t)| 4 , (210) 



However, to gain more physical insight into the solution of the Gross-Pitaevskii equation it 
is helpful to apply the collective coordinate method mentioned above. This is in particularly 
true because for a gas with effectively attractive interactions the breathing mode is the 
only unstable mode of the condensate ||115| , |119|| . Therefore, it appears to be an excellent 



approximation to the effective action Se[p] to expand up to quadratic terms in the coefficients 
Cn£ m (t), but to take all orders of q{t) into account. Clearly such an approximation is 
not possible if we describe the breathing mode with an additional coefficient Cioo(^) as in 

Eq. dog). 

To complete our discusion of the collisionless dynamics of a condensate with negative 
scattering length, we thus need to add to the noninteracting action in Eq. (|198|) the various 
contributions from the interaction term 

2nati 2 rW 



m 



f p f 

J drj dxp (x, r) 



in the action Se[p\- In lowest order in the coefficients C r ^ m {j) this only adds the mean-field 
or Hartree term — NQ?i 2 \a\/ \/2nmq 3 to V(q) as we have already discussed in detail above. 
The contributions linear in C n £ m (r) are 

4irah 2 ^„ fW fr f Xoo(x;g) 

— ? Jo dT Cn °° J dX P(X; q) W = 



4nah V" ^ A n 1 (-1)^(271 + 1)! 
°^T^ Jo °" 00 Nf 2^7^! 
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Note that we have only linear terms for the coefficients C n oo( r )- This is as expected because 
in the metastable state the average density profile of the condensate is given by 

(p(x)> = p(x; (q)) + J2"(C nem ) %^§ , (211) 



■nlrn 



7T 



( g )2)3/4 



where due to the rotational symmetry of the isotropic harmonic oscillator potential the 
average value of the coefficients C n £ m (r) must be zero if the angular momentum variables 
(£m) 7^ (00). The average value of C n oo(T) can be calculated if we know the quadratic terms 
in C n oo(T), which become 

2nah 2 „ 1 T(n + n' + 3/2) 

'2.^ t^nOO^n'OO: 



(V) 3/2 ,/27r(2n + l)!(2n / + l)! 



n,n' 

For the calculation of all the static properties of the condensate and in particular a more 
accurate value of the maximum number of atoms N m , the latter two contributions are all 
that we need. However, for a study of its dynamics we require the full quadratic expression 



n,n' Im ^ \ H ) 

2Trah*i l ^rW, _ „ 1 T(n + n' + £ + 3/2) 



m ^ (Trg 2 ) 3 / 2 2 n+n>+e+3/2 ^ n \ T (n + £ + 3/2)n'!r(n' + £ + 3/2) 

As an application of the above approach, we consider the decay rate of the condensate 
due to thermal fluctuations. This rate equals ||106|| 

r(iV ,T) = ^exp(-^l , (212) 



2tt 1 I k B T 

where AV^ represents the macroscopic energy barrier for the condensate. Calculating the 
average (C n £ m ) as a function of q and substituting the result back into the action Se[p], 
we obtain a more accurate potential V(q) that now incorporates the fact that the average 
density profile of the condensate is not gaussian near the point of instability. With this 
more accurate potential we can again calculate the frequency of the 'breathing' mode co>ooo 
and the macroscopic energy barrier A V as a function of the number of condensate particles. 
For the experimental conditions of interest the results are given in Fiq. ^. We see that the 
maximum number of condensate particles is N m ~ 1316 in this case, which corresponds to 
N m ~ 0.60//|a| and should be compared to the exact result N m ~ 0.57//|a| obtained by 



Ruprecht et al. [104]. The difference is due to the fact that we have included in the action 
only terms up to quadratic order in C n £ m (r). By including higher order terms we should in 
principle be able to reproduce the exact value of N m , although this will be rather tedious in 
practice. The advantage of the present approach is that we can systematically incorporate 
the effect of fluctuations of the condensate, which cannot be achieved by a solution of the 
Gross-Pitaevskii equation. Hereby, we conclude our treatment of the dynamics of Bose- 
Einstein condensation in an interacting Bose gas with negative scattering length. From now 
on we will again focus on a gas with effectively repulsive interactions, in which case the 
experiments are in general in the more complicated strong-coupling limit. 
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C. Phase 'Diffusion' 



As we have just seen explicitly, a particularly interesting consequence of the finite size 
of the gas is that quantum fluctuations play a much more important role. Although this is 
especially true for the case of attractive interactions, it is also true for a Bose gas with repul- 
sive interactions. A striking example in this respect is the phenomenon of phase 'diffusion', 
which was recently discussed by Lewenstein and You [|136|| . We rederive their results for a 
trapped Bose gas in a moment, but first consider also the same phenomenon for a neutral 
and homogeneous superconductor. Not only is it possible in this manner to bring out the 
physics involved more clearly, but we are then also able to point out an important distinction 
between a fermionic and a bosonic superfluid that is also crucial for a good understanding 



of the concept of a quasicondensate in weakly-interacting Bose gases |24 . 

At zero temperature the dynamics of the superconducting order parameter, i.e., the BCS 
gap parameter A(x, t) that is proportional to the wavefunction of the condensate of Cooper 
pairs, is in a good approximation determined by a time-dependent Ginzburg-Landau theory 



137 13£ 1 with the action 



S[A*,A] 



N(0) 



dt I dx. 
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(213) 



where N(0) is the density of states for one spin projection at the Fermi energy ep = mv F /2 
and Ao is the equilibrium value of the order parameter [|140|| . Writing the complex order 
parameter in terms of an amplitude and a phase, we immediately observe that the ampli- 
tude fluctuations are gapped ||141|1 and can, therefore, be safely neglected at large length 
scales. The long-wavelength dynamics of the superconductor is thus dominated by the phase 
fluctuations, according to the action 




S[9) 



This also implies that the global phase 6o{t) 
dynamics that is governed by 



3 



(vef 



(214) 



/ (ix 6*(x, t)/V of the superconductor has a 



(215) 



using the fact that the total volume V of the system is given by N/n. 

Up to now our discussion has again been semiclassical. To consider also the quantum 
fluctuations, we have to quantize this theory by applying the usual rules of quantum me- 
chanics. Doing so, we find that the wavefunction of the overall phase obeys a Schrodinger 
equation 



n 



o 2 



N(0)N 891 



(216) 



with a 'diffusion' constant that can easily be shown to be equal to (2/h)dep/dN [p.42|| and 
is, most importantly for our purposes, proportional to 1/N. In the thermodynamic limit 
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iV — > oo a state with a well defined stationary phase is clearly a solution and we are then 
dealing with a system having a spontaneously broken £7(1) symmetry. However, for a finite 
(and fixed) number of particles the global phase cannot be well defined at all times and 
always has to 'diffuse' in accordance with the above Schrodinger equation. Note also that in 
the groundstate the phase is fully undetermined and |^(^o; ^)| 2 = l/2n. Maybe surprisingly, 
the same calculation is somewhat more complicated for a Bose gas because the amplitude 
fluctuations of the order parameter cannot be neglected even at the largest length scales. 
However, taking these amplitude fluctuations into account properly, we nevertheless arrive 
at an action that is equivalent to Eq. ( |215[ ) and hence again leads to the phenomenon of 
phase 'diffusion'. 

We start again from the action Se[p, d; p] for the condensate. The difference with the 
previous subsection is, however, that now we are not so much interested in the dynamics of 
the density but in the phase dynamics instead. Therefore we now want to integrate over the 
density field p(x, r). This cannot be done exactly and we therefore here consider only the 



strong-coupling limit, which was also treated by Lewenstein and You |[L36|1 . In that limit 
we are allowed to neglect the gradient of the average density profile [ 107|| and the action 
Se[p, 0; fj] is for the longest wavelengths well approximated by 



S E [p,0;p] = dr J dx ^p(x, T )—^ 



00(x, r) 

2-xah 2 



+V^(x)p(x, r) - pp(x, r) + p^(x, r) . (217) 

m J 

In equilibrium the average density profile of the condensate thus obeys 

<P(x)> = (A* " V e %*)) 6(p - V~(x)) . (218) 

Performing now the shift p(x, r) = (p(x)) + 5p(pc, r), we find for the zero-momentum part 



of the action 143 



S E [8N , 6 ; p] = hpEofa) + £ P dr {iUN^ + ^^ji^o) 2 ) , (219) 
where E (pi) and V (p) correspond, respectively, to the energy and the volume of the conden- 



sate in the so-called Thomas- Fermi approximation 1 144|| . Moreover, 5N (t) = J <ix <5p(x, r) 
represents the fluctuations in the total number of condensate particles in that same approx- 
imation, implying that the density fluctuations <5p(x, r) are only nonzero in that region of 
space where the condensate density does not vanish. 

Performing now the integration over the number fluctuations 5N (t) and the usual Wick 
rotation to real times, we immediately see that the effective action for the global phase of 
the condensate has precisely the same form as in Eq. ( |215|) , i.e., 



mVo(p) f ( d0 t 



2 



s ^^f;*(^J • (220) 

The appropriate 'diffusion' constant is therefore equal to 2irah/mVo(fi), which can easily be 
shown to be equal to (l/2h)dp/ 8Nq if we make use of the fact that in the Thomas- Fermi 
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approximation the chemical potential obeys /i = mw 2 R\- F /2 and the radius of the condensate 
is given by Rtf = (lSa^A^o/m 2 ^ 2 ) 1 / 5 [[L28|| . Hence, the 'diffusion' constant is proportional 
to 1/N . Note that if the condensate where contained in a box the 'diffusion' constant 
would be proportional to l/N instead. It is important to note also that, in contrast to the 
case of a fermionic superfluid, we have to integrate over the amplitude fluctuations of the 
order parameter to arrive at a quadratic action for the phase fluctuations. This leads to the 
important conclusion that for a bosonic superfluid it is impossible to be in a state with only 
phase fluctuations and no density fluctuations, even at the largest length scales. A more 
implicit way to arrive at the same conclusion was in fact already presented in Sec. |V[ since 
it is precisely the coupling between the amplitude and phase fluctuations that is required to 
understand that the dynamical exponent z is equal to 1. 



D. Strong-Coupling Limit 

The difficulty with the strong-coupling limit is that the appropriate one-particle states 
of the gas are no longer the states x«( x ) of the trapping potential because these states 
are now strongly coupled by the mean-field interaction. For our Caldeira-Leggett model 
this problem could be easily resolved by introducing the states Xa( x ) that incorporate the 
average interaction with the reservoir. However, for a trapped Bose gas the same procedure 
would lead to an expansion into states x' Q (x;t) that parametrically depend on time since 
the mean-field interaction depends on time also. In principle, we thus have to deal with the 
situation that not only the occupation numbers vary with time but also the states to which 
these occupation numbers belong. This results in highly nontrivial dynamics since it is a 
priori clear that an adiabatic approximation, in which we neglect the nondiagonal matrix 
elements / dx x'a ( x > ^)(^/^)xL'( x ! 0> ^ s n °t van d, because of the lack of a separation of 
time scales. Notice that we were able to circumvent the above problems for a Bose gas with 
effectively attractive interactions due to the fact that on the one hand the growth of the 
condensate essentially takes place in the weak-coupling regime and that on the other hand 
the collapse, which certainly occurs in the strong-coupling regime, requires in first instance 
only knowledge of the collisionless dynamics of the condensate. 

Notwithstanding the above remarks, we can make progress in the strong-coupling limit 
if we realize that in this limit the typical length scale on which the (noncondensate) density 
of the gas varies is always much larger than the trap size. We are therefore justified in 
performing a gradient expansion. To do so in a systematic fashion, it is convenient to 
express the various selfenergies and interactions in the coordinate representation, i.e., by 
introducing 

^W^(x,x';e) = ^X a (x)fisS K (e))(:,(x') (221) 

a,a' 

and similarly 

VW(x,x';y,y' ;e ) = £ E X^xMV^^xlWxW) ■ (222) 

a,a' P,/3' 

In this coordinate representation, our Fokker-Planck equation acquires in all generality the 
form 
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dt 
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50(x)50*(x') 



+ V e *(x) - 5(x - x') + SE< + »(x, x')) 0(x) 

4- / dydy' 7< + »(x, x'; y, y')0*(y)0(y W) j P[<f , 0; t] 
+ V e:c (x) - 5(x - x') + ftE^x',*)) 0*(x') 
+ J dydy' \/(-)(x',x;y',y)0*(x')0*(y')0(y) 
ft£ K (x, x') 

+ J dydy' ^(x,x';y,y')0*(y)0(y')} J P[0*,0;t] - (223) 



To study its physical content and to familiarize ourselves with how to perform the desired 
gradient expansion, we first analyze again the normal state of the gas in which we can neglect 
the nonlinear terms representing the effect of the condensate. 

In the normal state the condensate is absent and we are only interested in the behavior 
of the quantity (0(x)0*(x'))(t). From the above Fokker-Planck equation we see that it obeys 
the equation of motion 



^<0(x)0V)>W 



dx" ( 
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( h 2 d 2 



2m d. 



+ V™(x) - Ai(t)j 5(x - x") + ftE<+)(x, x") j (</>(x>* (x'))(t) 
^ + V-(x') - - x') + fiEH^^j (0(x)0*(x")>(t) 



--fcE^fx.x') . 
2 v ' 1 



Introducing now the Wigner distribution ||145|| by means of 



iV(x, k; t) + - = / dx' e~ lkx (0(x + x'/2)0*(x - x'/2))(t) 



(224) 



(225) 



and similarly the selfenergies hH^' K (x, k) by 

KW'^x, k) = /" dx' e- ik ' x '^S( ± ^(x + x'/2, x - x'/2) , (226) 
we find after substituting these relations in lowest nonvanishing order in the gradient d/dx 



the expected quantum Boltzmann equation [p40| , |53 



d Ar/ , , lde'(x,k;t) a . , , I9e'(x,k;t) 9 Ar/ , , 

-T out (x, k; £)iV(x, k; t) + T in (x, k; t)(l + N{x, k; t)) , (227) 
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where the energies of the atoms have to be determined selfconsistently from 



e'(x, k; t) = e(k) + V ex {x) + Re[hZ w (x, k; e'(x, k; t) - //(*))] . (228) 

Moreover, the expressions for the retarded energy ftE(+)(x,k;e'(x,k;f) - //(£)), and the 
collision rates r ou '(x, k; t) and r m (x, k;t) turn out to be identical to the ones derived for 



a homogeneous gas in Sees. [1 1.1] and |IV] We only need to replace e'(k; t) by e'(x, k; £) and 
N(k; t) by iV(x, k; t). Notice also that the same substitution has to be made in the Bethe- 
Salpeter equation for the many-body T matrix, which implies that the effective interaction 
k") now implicitly depends both on the position x and the time t. 
As an example, let us consider the gas not too close to the critical temperature. Then 
the effective interaction is well approximated by the two-body T matrix and we simply have 

e' x, k; t) = e k + V ex x + . 229 

m 

As a result the equilibrium solution iV(x, k; oo) of the quantum Boltzmann equation equals 
a Bose distribution evaluated at e(k) + V ex (x) + 87m(x)ah jm — \i. The equilibrium density 
profile n(x) = J dk N(pc, k; oo)/(27r) 3 thus corresponds to a local-density approximation but 
is indeed very accurate for the atomic gases of interest. Apart from describing the relaxation 
of the gas towards equilibrium, the quantum Boltzmann equation can also be used to study 
the collective modes of the trapped gas above the critical temperature. Using the same 
methods that have recently been applied to inhomogeneous Bose gases by Zaremba, Griffin, 
and Nikuni [ 146|| , we can for example derive the hydrodynamic equations of motion the gas. 



Moreover, the collisionless modes are obtained from 

—N x, k; t + — ■ j-N(x, k; t - -— V ex {x) + K -LJ • — N x, k; t ) = . 

at m ax hox \ m J ok 

(230) 

Solving this equation is equivalent to the random-phase approximation ||53|1 . More important 
is that the above collisionless Boltzmann equation can be shown to satisfy the Kohn theorem 
147|| by substituting the shifted equilibrium profile iV(x — x cm (t), k— (rn/fi)dx. cm (t)/dt] oo). 



For an harmonic external potential the center-of-mass of the gas cloud then precisely follows 
Newton's law <i 2 x cm (t) / 'dt 2 = — W ex (~x cm (t)) , which shows that the gas has three collision- 
less modes with frequencies that are equal to the three frequencies of the external trapping 
potential. 

The full quantum Boltzmann equation in Eq. Q227|) also describes the relaxation of the gas 
towards equilibrium. It can thus be used to study forced evaporative cooling by adding an 
appropriate loss term to account for the escaping high-energetic atoms. In agreement with 
our previous picture we then expect that, if the evaporative cooling is performed slowly, 
the distribution function iV(x, k; t) can be well approximated by a truncated equilibrium 
distribution with a time-dependent temperature and chemical potential. As a result, we 



can use the results from Sec. [TTT] to show that in the critical regime the gas again becomes 
unstable towards Bose-Einstein condensation. This is the case when e'(x, 0;i) < /i(t) for 
positions near the center of the trap where the density of the gas is highest. If that happens 
we are no longer allowed to neglect the nonlinear terms in the Fokker-Planck equation and 
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some corrections are neccesary. Applying again the Hartee-Fock approximation, which is 
correct sufficiently close to the critical temperature where k B T ^> n (x, t)T^ + \0, 0; 0), we 
find first of all that the quantum Boltzmann equation for the noncondensed part of the 
gas has the same form as in Eq. fl227|) . We only have to add the mean-field interaction 
2n (x, t)T^ + \0, 0, 0; 0) due to the condensate to the renormalized energy e'(x, k; t) and in 
the same way as in the homogeneous case also account for the presence of the condensate 
in the scattering rates r ou *(x, k; t) and T m (x, k; t). 

In addition, the resulting quantum Boltzmann equation for the noncondensed atoms is 
coupled to a Fokker-Planck equation for the condensate. After a gradient expansion the 
latter equals 

d 

ih—P o [(f)* o ,0 o ;t] = 

-/ <ix ^)(-l^ ! + l '" (x) -^ + fiS(+, ( x ' > 

+ T (+ »(0,0,0; 0)|^(x)| 2 j 0o(x)PoK, <k; t] 

+ / dx ^x)(- ! £ !+ ^ (x) -' iw+fiEM(x - o) 

+ T (+ »(0,0,0; 0)|A,(x)| 2 Wx)P„[,ft fo; i] 

1 

2 J b(po{x)b(po{x) 

where fiEW^fx.O) denotes hE&> K fak = 0). In more detail we have for the real part 
of the retarded and advanced selfenergies that 5*(x, 0) = Re[fi,E^(x, 0; e"(x, 0; t) — fi(t))]. 
Moreover, using our picture of the evaporative cooling of the gas, the imaginary parts follow 
most conveniently from the fluctuation-dissipation theorem 

*i?(x, 0) = -^W (X; 0) {-^~ + e'(x, 0; t) - p(t) + T< + >(0, 0, 0; O)|0 o (x)| 2 ) (232) 

where hT, K (x, 0) = hH^ ,K (x, 0; e"(x, 0; £)—//(£)) can in a good approximation be considered 
as independent of |0 O ( X )|. Near equilibrium the Fokker-Planck equation thus simplifies to 

2V72 



+ rW(0, 0, 0; O)|0 o (x)| 2 ) M*)Pol\M,t} 
4/dx^(x,0)^ y (-^ + e '(x ) 0)-, 

+ T«(0, 0, 0; O)|0 o (x)| 2 l 05(x)P o [|0o|; *] 
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-\ I dKKE K (i,i,0) — 6 ] P [\M,t] (233) 

and the probability distribution of the condensate relaxes to the expected stationary state 
P o [|0o|; oo] oc 

exp Lp J rfx 0*(x) (-^ + e'(x, 0) - p + gi °) | 0o(x )|A 0o(x )| . 



From this result we conclude that Eqs. ( f227| ) and ( |231| ) give an accurate discription of the 
dynamics of Bose-Einstein condensation in the strong-coupling regime. Unfortunately, the 
solution of these coupled equations is not possible analytically and we therefore have to 
refer to a future publication for a numerical study and a comparison of this theory with 
experiment. We here only remark that the semi classical approximation in this case amounts 
to the Thomas-Fermi result for the condensate density 

n ° (x; t] = rw^o o-o) " e ' (x ' 0; t]) 0(Mt) " e ' (x ' 0; t)} • (234) 

Such an approximation implicitly assumes, however, that the dominant mechanism for the 
formation of the condensate is by means of a coherent evolution and not by incoherent 
collisions. This might be expected not to be very accurate on the basis of the homogeneous 
result that the second mechanism only diverges logarithmically with the system size, whereas 
the first diverges linearly. Nevertheless it seems that only a full solution of Eqs. (|227|) and 
( |231|) can settle this important issue. 

E. Collective Modes 

Although we have up to now primarily been focussed on the dynamics of Bose-Einstein 
condensation itself, it is interesting that we have in fact also arrived at an accurate discription 
of the collective modes of a condensed gas at nonzero temperatures. Since these collective 
modes are receiving much attention at present, we now briefly want to discuss how such 
a description emerges from our Fokker-Planck equation |[148|| . Since most experiments at 
present are performed in the collisionless limit, we consider only this limit here. It is in 



principle possible, however, to derive the hydrodynamic equations of motions as well ||146 
Let us for simplicitly first assume that the temperature of the gas is sufficiently far below the 
critical temperature that the many-body T matrix is well approximated by the two-body 
T matrix, but not so low that the Hartree-Fock approximation is no longer valid. In this 
regime we have that ^S^(x, k) = 87m'(x, t)ah 2 /m and the dynamics of the noncondensed 
part of the gas is again determined by the collisionless quantum Boltzmann equation given 
in Eq. ( |230| ), where the total density now consists of two contributions 

n(x, t) = no(x, t) + n'(x, t) = | (0(x)) (t) | 2 + / -— - iV(x, k; t) . (235) 



(2tt) 3 

The dynamics of the condensate follows again from the nonlinear Schrodinger equation 
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d f h 2 V 2 Airah 2 1 

? ^_(0(x))(t) = |-^ + \/-(x) -/,+ —— (2n'(x,t) + n o (x,t))j(0(x)>(t) , (236) 

which can be rewritten as the continuity equation 

-|n (x, t) = -V ■ (v s (x, £)n (x, £)) (237) 
together with the Josephson relation for the superfluid velocity 

JU s (x, t) = -V (v ex {x) ^-(2n'(x, t) + no(x, £)) + imv s 2 (x, ^ . (238) 



Here, we used that (0(x))(t) = yno(x,t)e 4 ^ x,f ', v s (x, t) = W6 l (x, t)/m, and have applied 
the Thomas- Fermi approximation, which is appropriate in the strong-coupling limit. It will 
presumably not come as a surprise that these coupled equations for iV(x, k; t), n (-x, t), and 
v s (x,i) can be shown to fulfill the Kohn theorem exactly. This is an important constraint 
on any theory for the collective modes of a trapped gas, as stressed recently by Zaremba, 
Griffin, and Nikuni 111461 • 



As we have mentioned previously, in trapped atomic gases it is possible to cool the gas 
even to such low temperatures that k B T <C n (x; t)T^ +S) (0, 0; 0). In this regime it is no longer 
allowed to use the Hartree-Fock approximation and we must use the Bogoliubov approxi- 
mation instead. This is achieved as follows. In the Hartree-Fock approximation we have 
taken for the probability distribution P[(fi*, <fi;t) the separable form Po[0o> 'Po', t)Pi[(f>'*, <f)'; t). 
For such an ansatz the U(l) invariance of the theory requires in fact that P[0*,0;i] is 
invariant under independent phase transformations on <f) and <p'. Clearly this ansatz is 
therefore too symmetric in general, since it is only required that P[<f>*,<f>',t] is invariant 
under simultaneous phase transformations on <f) and <$' . For such more general solutions 
the probability distributions for the condensed and noncondensed parts of the gas are then 
given by P a [^M = /<#'*]<#'] PWA'A and Pi[0'*,0';t] = M0»q] P{^,<P;t], re- 
spectively. The Fokker-Planck equations for these probability distributions can thus be 
found from integrating the general Fokker-Planck equation in Eq. ([223 ) first over 0'(x) and 



subsequently over 0o( x )- Taking care that we do not double count effects of the interaction 
that are already accounted for in the various selfenergies, we find after a gradient expansion 
that the Fokker-Planck equation for the condensate is identical to Eq. ( [23 1| ). The collision- 
less dynamics of the condensate is therefore again determined by the nonlinear Schrodinger 
equation in Eq. ( |236| ) or in the Thomas- Fermi approximation by Eqs. ( |237] ) and ( |238| ). The 



Hartree-Fock and Bogoliubov approximations are clearly identical in this respect. This is, 
however, not true if we consider the dynamics of the noncondensed part of the gas. 
Integrating over 0o( x ) we now obtain in the collisionless limit that 

" / dx WW { (~ST + v (x) - " + — n(x ' () j * (x) 

1 ;; " fi ^„(x, f )e 2 »^V>)l WW; «] 



rn 
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+ / dx W¥) { r^r + v (x) - " + — n(x - 4) J * (x) 

+^^n„(x,t)e-»W>0'(x)|p 1 [r>';«] ■ (239) 

This implies that, in contrast to the Hartree-Fock treatment, the equation of motion for 
(0'(x)0'*(x'))(t) is coupled to the equation of motion of the so-called anomalous average 
(0'*(x)0'*(x / ))(t). In the framework of a qradient expansion, these equations can however 
be decoupled by performing a local Bogoliubov transformation |17|]. In this manner we can 
then show that the Wigner distribution for the Bogoliubov quasiparticles obeys the expected 
collisionless Boltzmann equation 

m N(x, k; t ) + -JgL-i • ^(x, k; t) ^ . -7V(x, k; t ) = , (240) 



with the dispersion given by 



Snah 2 



kuj(x, k; t) = \/e 2 (k) + n (x, t)e(k) + hk ■ v s . (241) 



ID 



Note that this shows that the distribution function N(pc, k; t) gives the number of Bogoliubov 
quasiparticles with momentum hk at position x in the frame in which the superfTuid velocity 
vanishes |14],^|. Note also that the above derivation of the Bogoliubov theory is particle- 
number conserving since the probability distribution P[<f>*, 4>; t] is explicitly U(l) invariant. 
It is therefore equivalent to the approach recently put forward by Gardiner [ |149|| , and Castin 



and Dum ||150|| . More important for our purposes is that the coupled equations for iV(x, k; t), 
no(x, t), and v s (x, t) can again be shown to fulfill the Kohn theorem exactly. We therefore 
believe that they give an accurate description of the collisionless modes of a trapped Bose 
condensed gas, that may be able to explain the existing discrepancy between theory and 
experiment. 

As a first step towards this goal we here consider the time-dependent Hartree-Fock theory 
for the collective modes of the Bose condensed gas, in which the collisionless Boltzmann 



equation in Eq. (|230|) is coupled to the nonlinear Schrodinger equation in Eq. (p36|) . It 



turns out that a numerical determination of the eigenmodes of the collisionless Boltzmann 
equation is rather difficult. Therefore, we have recently proposed a variational approach to 
solve these coupled equations [ |151| . Denoting the anisotropic harmonic trapping potential 



used in the experiments ||123| , |124|| by V ex (x) = J2j muJ j x 'j/2, we first of all take for the 
condensate wavefunction the gaussian ansatz 

with the functions qj(t) describing the 'breathing' of the condensate in three independent 
directions. Clearly, this is the anisotropic generalization of the gaussian ansatz that we 
used in Sec. |V B| to study the stability of a condensate with negative scattering length. 
Although it is not immediately obvious that such an ansatz is also appropriate in the strong- 



coupling regime of interest here, it has been shown by Perez-Garcia et al. |118|| that at zero 
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temperature the assumption of a gaussian density profile nevertheless leads to rather accurate 
results for the frequencies of the collective modes. We therefore anticipate that at nonzero 
temperatures the same will be the case. 

Next we also need an ansatz for the distribution function iV(x, k; t) of the thermal cloud. 
In Ref. 



151] we considered the gaussian 



iV(x,k;t) = iV'n 



X 



ri ex p 




da.j{t) 



ha>j(t) dt 



+ 



3 



2 aj(t) 



(243) 



where otj(t) are again three scaling parameters, N' = N — N is the number of noncondensed 
particles, and (otj) denotes the equilibrium value of otj{t) that in general is greater than 1 
due to the repulsive interactions between the atoms. Here, we however also consider the 
more accurate Bose distribution 



iv( X ,M=iv'-— n 



flUlj 



C(3) 



k B T{a j ) 2 ) 
~h 2 a 2 At) 



k B T(aj) 2 



2m 



kj 



mxj 



dotj {t) 



%a.j(t) dt 



+ 



X ; 



2 a](t) 



(244) 



with ((3) ~ 1.202. In both cases the functional form of the distribution function is physically 
motivated by the following reasons. First, by integrating over the momentum Kk we see that 
the time-dependent density profile of the noncondensed cloud n'(x, t) is a time-independent 
function of Xj/otj(t), as desired. Second, for such a density profile we must have that the 
local velocity h{kj){x,t)/m is given by (xj/otj(t))dotj(t)/dt. Third, we want the equilibrium 
distribution to be only anisotropic in coordinate space, and not in momentum space. This 
explains the appearance of the factor 1/ (a-,) 2 in the exponent. 

Using the ideal gas result for the number of condensate particles and substituting the 
above ansatz for the condensate wavefunction and the Wigner distribution function into the 
nonlinear Schrodinger and collisionless Boltzmann equations, it is straightforward but some- 
what tedious to derive the six coupled equations of motion for the 'breathing' parameters 
qj(t) and atj(t). Linearizing these equations around the equilibrium solutions (qj) and {otj), 
respectively, we can then also obtain the desired eigenmodes and eigenfrequencies. The re- 
sults of this calculation for the trap parameters of Jin et al. ||123|| are shown in Fig. |5] together 



with the experimental data. Note that the trap has a disk geometry with uj x = uo y = uo r 



and uj z = y%uj r . Due to the axial symmetry, the modes of the gas can be labeled by the 
usual quantum number m that describes the angular dependence around the z-axis. On 
the whole there appears to be a quite reasonable agreement between this simple theory and 
experiment. The only exceptions are the two data point almost halfway in between the in 
and out-of-phase m = modes. A possible explanation for this discrepancy is that exper- 
imentally both modes are excited simultaneously ||152|] . To make sure of this, however, we 



need to determine the oscillator strength of the various modes. Work in this direction is in 
progress and will be presented in a future publication. 
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VI. DISCUSSION AND CONCLUSIONS 



The purpose of this paper has been twofold. First, we have derived a single Fokker- 
Planck equation for the nonequilibrium dynamics of a weakly interacting Bose gas. This 
Fokker-Planck equation is capable of describing simultaneously the coherent and incoherent 
processes taking place in the gas and can be used for homogeneous as well as inhomogeneous 
situations. It can be derived in a very natural way in the Schwinger-Keldysh formalism, to 
which we have therefore also tried to give an introduction. The Fokker-Planck equation is 
here primarily used to study the dynamics of Bose-Einstein condensation in homogeneous 
and inhomogeneous atomic Bose gases. In particular, by treating both the weak-coupling 
and strong-coupling limits of the theory, we have been able to consider not only trapped 
atomic Bose gases with effectively attractive but also with effectively repulsive interactions, 
which leads to a better understanding of the fundamental difference between the physics 
in these two cases. Although we have only applied the Schwinger-Keldysh formalism to 
atomic Bose gases here, it can also be used to study nonequilibrium problems in trapped 
Fermi gases. An example is the study of the timescale for the formation of Cooper pairs in 
atomic 6 Li ||153|| , which due to the anomalously large and negative value of the interatomic 
scattering length for this species, is of considerable interest at present. 

Second, and maybe most importantly, with our Fokker-Planck equation we have been 
able to make some theoretical progress on three topics in the field of Bose condensed atomic 
gases where there exists at present a discrepancy between theory and experiment. We 
refer here to the experiments on the formation of the condensate, on the spectrum of the 
collisionless modes and on the collapse of a condensate with negative scattering length. The 
main difficulty with the first topic is that the experiments are performed in the strong- 
coupling limit. In the language of quantum optics this implies that we are dealing with 
a multi-mode situation. Although this complicates the theory considerably, it also makes 
the physics involved much richer because coherent interaction effects now compete with the 
incoherent collisions. In this regime we may therefore hope to study in some detail the 
formation of topological defects in the wake of a phase transition, which turns out to be an 
important problem in cosmology as well |154j|. Moreover, recent experiments have shown 



that by a sinusoidal modulation of the trapping potential it is possible to adiabatically 
cycle through the phase transition | |155| |. A theoretical description of this experiment is 
easily obtained within our strong-coupling formalism. We only have to replace V ex (x) by 
V ex (x,t). 

The second discrepancy between theory and experiment is associated with the eigenfre- 
quencies of the collisionless modes of a Bose condensed gas at nonzero temperatures, when 
also a considerable density of noncondensed atoms is present in the trap. The most accu- 
rate theory that has recently been applied to this problem, is an approximate many-body 
T-matrix calculation in which the noncondensate part of the gas is assumed to have no 
dynamics ||156|| . Due to this last feature the modes of the gas are essentially the modes of 



a condensate in the effective potential V ex (x) + 2n'(x)T( + )(0, 0, 0; 0), which near the center 
of the trap has effectively smaller harmonic oscillator frequencies than the external trapping 
potential. Therefore, the Kohn theorem is not satisfied. As we have shown, this problem 
can be overcome if the dynamics of the noncondensed part of the gas is determined by an 
appropriate collisionless Boltzmann equation. Furthermore, by including the collision terms 
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we can also find the damping of the long-wavelength modes of interest in a similar way as 



Kavoulakis, Pethick, and Smith have done for the normal state of the gas ||157 



The last theoretical challenge that we have touched upon is an accurate determination 
of the number of atoms that remain after a single collapse of a condensate with nega- 
tive scattering length. We have argued that this in principle requires the solution of the 
Gross- P it aevskii equation in combination with a rate equation for the number of condensate 
particles. In addition, we have indicated how we can find an approximate description of the 
dynamics of the collapse that systematically improves upon our gaussian results. As in the 
two previous problems, work is in progress to numerically investigate if the theory presented 
in this paper is indeed sufficiently accurate to resolve the discrepancies with experiment. 
Although we cannot be sure about the outcome of these studies at present, we believe that 
the many-body T-matrix approximation used throughout this paper is an excellent starting 
point in all three cases, and that our Fokker-Planck approach to the dynamics of weakly in- 
teracting Bose gases presents a convenient framework in which to derive in a straightforward 
way also the necessary corrections if they turn out not to be negligible. 
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FIGURES 




FIG. 1. Semiclassical picture of the coherent stage of Bose-Einstein condensation 
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FIG. 2. Renormalization group calculation for the (meta) stable phases of a homogeneous Bose 
gas with two and three-body interactions. The two-body T matrix is equal to 47rh 2 a/m, whereas 
the three-body T matrix equals A-K?i 2 b A /m. 



76 




77 




FIG. 4. Dynamical properties of a Li condensate. Curve 1 gives the energy barrier AV for 
thermal fluctuations and curve 2 the frequency for the unstable breathing mode. Together these 
quantities determine the decay rate of the metastable condensate at the relatively high temperatures 
achieved in current experiments. 
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FIG. 5. Collsionless modes in a Bose condensed 8 7Rb gas. The dashed lines give the results 
with a gaussian ansatz for the distribution function of the thermal cloud and the full lines for a 
Bose distribution ansatz. Curves 1 and 2 correspond to the m = in and out-of-phase modes, 
respectively. Similarly, curves 3 and 4 give the frequencies of the m = 2 in and out-of-phase modes. 
The experimental data is also shown. Triangles are for a m = mode and circles for a m = 2 
mode. 
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